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Abstract 

We construct finite-dimensional invariant manifolds in the phase space of the 

■ Navier-Stokes equation on and show that these manifolds control the long-time 
. behavior of the solutions. This gives geometric insight into the existing results on 

^ I the asymptotics of such solutions and also allows one to extend those results in a 

■ number of ways. 

o ■ 

5 ■ 1 Introduction 

J3 In the last decade and a half, starting with the work of M. Schonbek, the long-time 

behavior of solutions of the Navier-Stokes equation (and the related vorticity equation) 
on unbounded spatial domains has been extensively studied. (See ||14|, |H1 H], 



^ , [|T^, ||2^, |21[| and pO|] for a small sampling of this literature.) This prior work used 

^ I a variety of techniques including energy estimates, the Fourier splitting method, and a 
detailed analysis of the semigroup of the linear part of the equation. In the present paper 
we introduce another approach, based on ideas from the theory of dynamical systems, 
to compute these asymptotics. We prove that there exist finite-dimensional invariant 
manifolds in the phase space of these equations, and that all solutions in a neighborhood 
of the origin approach one of these manifolds with a rate that can be easily computed. 
Thus, computing the asymptotics of solutions up to that order is reduced to the simpler 
task of determining the asymptotics of the system of ordinary differential equations that 
result when the original partial differential equation is restricted to the invariant manifold. 
Although it is technically quite different from our work, Foias and Saut, have also used 
invariant manifold theory to study the long-time behavior of solutions of the Navier-Stokes 
equation in a bounded domain. 

As we shall see, if one specifies some given order of decay in advance, say 0{t~''), one 
can find a (finite-dimensional) invariant manifold such that all small solutions approach 
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the manifold with at least that rate. Therefore, one can in principle compute the asymp- 



totics of solutions to any order with this method. In Subsection [4.1| , for instance, we 
illustrate how these ideas can be used to extend known results about the stability of the 
Oseen vortex. In particular, we show that as solutions approach the vortex solution their 
velocity fields have a universal profile. 

If one adopts a dynamical systems point of view towards the question of the long-time 
behavior of solutions of the Navier-Stokes equation, many phenomena which have been 
investigated in finite-dimensional dynamics immediately suggest possible effects in the 
Navier-Stokes equation. As an example, one knows that in ordinary differential equations 
the presence of resonances between the eigenvalues of the linearized equations may produce 
dramatic changes in the asymptotic behavior of solutions. Exploiting this observation in 



Subsection we construct solutions of the Navier-Stokes and vorticity equations whose 
asymptotic expansion contains logarithmic terms in time. 

Finally, we feel that the geometric insights that the invariant manifold method pro- 
vides are very valuable. As an example of the sorts of insights this method provides, we 
reexamine the results of Miyakawa and Schonbek |^ on optimal decay rates in Subsec- 



tion We find that the moment conditions derived in have a very simple geometric 
interpretation - they are the analytic expression of the requirement that the solutions lie 
on certain invariant manifolds. 

To explain our results and methods in somewhat more detail, we recall that the Navier- 
Stokes equation in has the form 

^ + (u ■ V)u = Au - Vp , V ■ u = , (1) 
ot 

where u = u(x, t) G is the velocity field, p = p{x, t) G R is the pressure field, and 
X G R^, t > 0. For simplicity, the kinematic viscosity has been rescaled to 1 in Eq.(|I]). 

Our results are based on two ideas. The first observation is that it is easier to derive 
the asymptotics of (|I]) by working with the vorticity formulation of the problem. That is, 
we set uj = rot u = diU2 — d2Ui and study the equation (^) for oj rather than (p itself. One 
can then recover the solution of the Navier-Stokes equation via the Biot-Savart law (^) 
(see Lemma and Appendix ^ for estimates of various norms of the velocity field in 



terms of the vorticity.) While there have been some studies of the asymptotics of solutions 
of the vorticity equation (notably 0] and the relationship between the vorticity and 



velocity does not seem to have been systematically exploited to study the asymptotics of 
the Navier-Stokes equations. We feel that approaching the problem through the vorticity 
equation has a significant advantage. It has been realized for some time that the decay 
rate in time of solutions of (|l]) is affected by the decay rate in space of the initial data. For 
instance, general solutions of the Navier-Stokes equation in L^(R^) with initial data in 
L^nL^ have L?' norm decaying like Ct"^/^ as t ^ +oo. However, it follows from Wiegner's 
result 1^ that, if the initial data satisfy /j^2(l + |2^|)|uo(a;)| dx < cxd, this decay rate can be 
improved to Ct~^. However, the semiflow generated by the Navier-Stokes equation does 



not preserve such a condition. For example, in Subsection [4.3| , we prove that there exist 
solutions of (|T]) which satisfy /j^2(l + |3;|)|u(a;, t)| dx < oo when t = 0, but not for later 
times. The vorticity equation does not suffer from this drawback - as we demonstrate in 
Section ^, solutions of the vorticity equation in R^ which lie in weighted Sobolev spaces 
at time zero remain in those spaces for all time. A similar result holds for small solutions 
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of the vorticity equation in R^, see [0. As we will see in later sections, these decay 
conditions are crucial for determining the asymptotic behavior of the solutions, and the 
fact that the vorticity equation preserves them makes it natural to work in the vorticity 
formulation. It should be noted, however, that this approach requires the vorticity to 
decay sufficiently rapidly as |x| — > oo and this is not assumed in |28|, for instance. On the 
other hand, assuming that the vorticity decreases rapidly at infinity is very reasonable 
from a physical point of view. This is the case, for instance, if the initial data are created 
by stirring the fiuid with a (finite size) tool. Note that in contrast, even very localized 
stirring may not result in a velocity field that decays rapidly as x tends to infinity. As an 
example, in Subsection we discuss the Oseen vortex which is a solution of (|l|) whose 
vorticity is Gaussian, but whose velocity field is not even in L^. 

The second idea that helps to understand the long-time asymptotics of (|l]) is the 
introduction of scaling variables (see (|T^) and (0) below). It has, of course, been realized 
for a long-time that in studying the long-time asymptotics of parabolic equations it is 
natural to work with rescaled spatial variables. However, these variables do not seem to 
have been used very much in the context of the Navier-Stokes equation. (Though they are 
exploited in a slightly different context in (0.) For our work they offer a special advantage. 
If one linearizes the Navier-Stokes equation or vorticity equation around the zero solution, 
the resulting linear equation has continuous spectrum that extends all the way from minus 
infinity to the origin. If one wishes to construct finite-dimensional invariant manifolds in 
the phase space of these equations which control the asymptotics of solutions, it is not 
obvious how to do this even for the linearized equation - let alone for the full nonlinear 
problem. However, building on ideas of 1 27] we show that, when rewritten in terms of the 
scaling variables, the linearized operator has an infinite set of eigenvalues with explicitly 
computable eigenfunctions, and the continuous spectrum can be pushed arbitrarily far into 
the left half plane by choosing the weighted Sobolev spaces in which we work appropriately. 
(See Appendix ^ for more details.) We are then able to construct invariant manifolds 
tangent at the origin to the eigenspaces of the point spectrum of this linearized equation 
and exploit the ideas that have been developed in finite-dimensional dynamics to analyze 
the asymptotic behavior of solutions of these partial differential equations. 

We conclude this introduction with a short survey of the remainder of the paper. 
Although our ideas are applicable in all dimensions we consider here the case of fluids 
in two dimensions. Because the vorticity is a scalar in this case many calculations are 
simplifled, and we feel that the central ideas of our method can be better seen without 
being obscured by technical details. Thus, in Section ^ we begin by surveying the (well 
developed) existence and uniqueness theory for the two-dimensional vorticity equation. 
In Section ^ we introduce scaling variables and the weighted Sobolev spaces that we use 
in our analysis. We then show (see Theorem |3.5|) that in this formulation there exist 
families of flnite-dimensional invariant manifolds in the phase space of the problem, and 
that all solutions near the origin either lie in, or approach these manifolds at a computable 
rate. Furthermore, we obtain a geometrical characterization of solutions that approach 
the origin "faster than expected" , see Theorem |3.8| . In Section ^ we apply the invariant 
manifolds constructed in the previous section to derive the above mentioned results about 
the long-time behavior of the vorticity and Navier-Stokes equations. In a companion paper 
T3| , we obtain similar results for the small solutions of the Navier-Stokes equation in three 



dimensions by a slightly different method. Finally, there are three appendices which derive 
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a number of facts we need in the previous sections. In Appendix we study the spectrum 
of the hnear operator L appearing in the rescaled vorticity equation (|I^) , and we obtain 
sharp estimates on the semigroup it generates. In Appendix ^, we study in detail the 
relationship between the velocity field u and the associated vorticity lo. In particular, we 
show how the spatial decay of u is related to the moments of the vorticity lo. Appendix |y 
derives a technical estimate on the invariant manifold constructed in Subsection 



Notation: Throughout the paper we use boldface letters for vector-valued functions, 
such as u(a;,t). However, to avoid a proliferation of boldface symbols we use standard 
italic characters for points in R^, such ). In both cases, I ■ I denotes the 

Euclidean norm in R^. For any p G [1, cxd] we denote by the norm of a function in 
the Lebesgue space L^'(R^). If f e (Lp(R^))^, we set = | |/| |p. Weighted norms play 
an important role in this paper. We introduce the weight function h : R^ R defined 
by b{x) = (1 + |a;p)2. For any m > 0, we set = l^'^/h, and denote the resulting 

Hilbert space by L'^{m). If / e C0([0, T], Lp(R2)), we often write f{-,t) or simply f{t) to 
denote the map x f{x,t). Finally, we denote by C a generic positive constant, which 
may differ from place to place, even in the same chain of inequalities. 
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Navier-Stokes equation. The research of C.E.W. is supported in part by the NSF under 
grant DMS-9803164. 



2 The Cauchy problem for the vorticity equation 

In this section we describe existence and uniqueness results for solutions of the vorticity 
equation. As stressed in the introduction, our approach is to study in detail the behavior 
of solutions of the vorticity equation, and then to derive information about the solutions 
of the Navier-Stokes equation as a corollary. The results in this section are not new and 
are reproduced here for easy reference. 

In two dimensions, the vorticity equation is 

+ (u ■ V)iu = Aiu , (2) 

where u = u{x, t) E H, x = (xi, X2) G R^, t > 0. The velocity field u is defined in terms 
of the vorticity via the Biot-Savart law 

^u;{y)dy, xeR". (3) 

y\ 

Here and in the sequel, if x = (xi, X2) G R^, we denote x = (xi, ^2)'^ and X"*- = (— X2, Xi)'^. 
The following lemma collects useful estimates for the velocity u in terms of 00. 
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Lemma 2.1 Let u be the velocity field obtained from uo via the Biot-Savart law 

(a) Assume that I < p < 2 < q < oo and i = i- i. If uo e Lp{H?), then u G ^^(R^)^, 
and there exists C > such that 

|u|g < C\u\p . (4) 

(b) Assume that l<p<2<q< oo, and define a G (0, 1) by the relation | = ^ + 
Ifu G Lp(R2) n ^^(R^), then u G L°°(R2)2, and there exists C > stxc/i i/iai 

|u|oo < . (5) 

(c) 74ss'ume that l<p<oo. //coiG Lp(R^), i/ien Vu G L^'(R^)'^ anc? i/iere exists C > 

|Vu|p < C|tj|p . (6) 
In addition, div u = and rot u = (9iM2 — 52^1 = uJ. 

Proof: (a) follows from the Hardy-Littlewood-Sobolev inequality, see for instance Stein 
||26| , Chapter V, Theorem 1. To prove (b), we remark that, for all i? > 0, 

|u(a;)| < — / \uj{x - y)\ — Ay + — I \uj{x-y)\— Ay 
27r J\y\<R \y\ 27r J\y\>R \y\ 

1-2 1 

< C\uj\qR 1 + C\uj\p , 

Rp 

by Holder's inequality. Choosing R = {\u\p/\uj\q)^ , where /3 = j^j^ = 2/p~i ' obtain 
(^). Finally, holds because Vu is obtained from u via a singular integral kernel of 
Calderon-Zygmund type, see [^, Chapter II, Theorem 3. QED 



The next result shows that the Cauchy problem for (|^) is globally well-posed in the 
space L^(R^). 

Theorem 2.2 For all initial data ujq G L^(R?), equation (j^ has a unique global solution 
uo G C°([0,oo),L1(R2)) n C°((0,oo),L°°(R2)) such that cj(0) = uoq. Moreover, for all 
p G [1, +oo], there exists Cp > such that 

Kt)|,<^^, t>0. (7) 

Finally, the total mass of uj is preserved under the evolution: 



uj{x, t)dx = j ujq{x) da; , t > . 

Proof: The strategy of the proof is to rewrite as an integral equation (see (|ll]) below), 
and then to solve this equation using a fixed point argument in some appropriate function 
space. We refer to Ben-Artzi [|| and Brezis [@ for details. QED 

Remark 2.3 L}{R?) is not the largest space in which the Cauchy problem for equation 
can be solved. For instance, it is shown in /|73|/, / |77l/ that ^ has a global solution 
for any uoq G M.{^), the set of all finite measures on R^. However, uniqueness of this 
solution is known only if the atomic part of ujq is sufficiently small. Moreover, such a 
solution is smooth and belongs to L^(R^) for any t > 0. Thus, since we are interested 
in the long-time behavior of the solutions, there is no loss of generality in assuming that 
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If uj{t) is the solution of given by Theorem p.2| , it follows from Lemma |2]l| that 
the velocity field u(t) constructed from u(t) satisfies u G C°((0, oo), L''(R^)^) for all 
q G (2, oo], and that there exist constants Cg > such that 

Ht)u<^4^, t>o. (8) 

Moreover, one can show that u[t) is a solution of the integral equation associated to (|ip. 

In the fact that we cannot bound the norm of the velocity field is not a technical 
restriction. As we will see, even if uj{x) is smooth and rapidly decreasing, the velocity 
field u(x) given by (0) may not be in L^(R^)^. A typical example is the so-called Oseen 
vortex, which will be studied in Section below. In fact, it is not difficult to verify 
that, if u G L^(R^)^ and u = rotu G L^^R"^), then necessarily J^2^{x) dx = 0. In this 
situation, the decay estimates (^, (H) can be improved as follows: 

Theorem 2.4 Assume that uq G L^(R^)^ and that ujq = rotuo G L^(R^). Let uj{t) be 
the solution of given by Theorem \2.2^ . Then 

lim t^~p|cj(t)|p = , l<p<oo. (9) 

Ifu{t) is the velocity field obtained from uj(t) via the Biot-Savart law then 

lim t5-||u(t)|^ = , 2<g<cx). (10) 



Proof: The decay estimate is a particular case of Theorem 1.2 in [^]. However, The- 
orem can also be proved by the following simple argument. Since VL{t) is a solution of 
the Navier-Stokes equation in L^(R^)^, it is well-known that Vu G -Z^^((0, +oo), L'^(R?Y). 
In particular, \uj{t)Wdt < oo. On the other hand, since ^|o^(t)|2 = -2|Vu;(t)||, the 
function 1 1— > |ix;(t)|2 is non-increasing. Therefore, 

t\uj{t)\l<2 f \uj{s)\lds'^= e{tf > 0, 

^ 2 

which proves (j^) for p = 2. 

We now use the integral equation satisfied by uj(t), namely 



Jo 

Since ujq G L-'^(R^) and Jj^a dx = 0, a direct calculation shows that |Lo'i(t)|i ^ as 
t — > +00. On the other hand, 

\u2{t)\i < I |u(s)u;(g)|idg < / |u(g)|2|a;(g)|2dg 

Jo Vt- s Jo y/t - s 



Vt- s a/s Jo a/x(1 -x) *^+°° 

by Lebesgue's dominated convergence theorem. Thus, (^ holds for p = 1, hence for 
1 < p < 2 by interpolation. Next, using the integral equation again and proceeding as 
above, it is straighforward to verify that (|^) holds for all p G [1, +oo]. From Lemma |2.1| , 
we then obtain (p!0|) for 2 < q < oo. Finally, the fact that |u(i(:)|2 — as t — +oo is 
estabhshed in [|l|. QED 
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3 Scaling Variables and Invariant Manifolds 



Our analysis of the long-time asymptotics of (§) depends on rewriting the equation in 
terms of "scaling variables" or "similarity variables" : 

If uj{x,t) is a solution of @ and if u{t) is the corresponding velocity field, we introduce 
new functions w{^,t), v(^,r) by 

^(^■*' = TT7'"(75^''°s(i + *)) ■ (12) 

Then w{^,t) satisfies the equation 

drW = Cw-{\- V^)w , (14) 

where 

C.tii = Af?/; -I- 

2 

and 



Cw = A^w + h^-V^)w + w (15) 



v(e,r) = ^ / ^-^^wiv,r)dv. (16) 

Scaling variables have been previously used to study the evolution of the vorticity by Giga 
and Kambe |14] and Carpio [^. 

The results of the previous section already provide us with some information about 
solutions of (|14D- For example, for all Wq G L^(R^), there exists a unique solution w G 
C0([0,oo),Li(R2)) n C0((0,oo),L°^(R2)) of such that w{0) = wq. Translating (0) 
and (H) into rescaled variables, we see that, for all r > 0, 

Hr)\, < ^y^, l<p<oo, (17) 
a^ry p' 

|v(r)|, < ^-^'^'l. , 2<g<oo, (18) 

where a(r) = 1 - e"^. Moreover, if vq = v(-, 0) G L^(R^)^, then /j^a ^^(C; t) d^ = for all 
r > 0, and it follows from (|) that |ti'(r)|p ^ as r ^ +oo for all p G [1, +oo]. 

However, in order to derive the long-time asymptotics of these solutions we need to 
work not only in "ordinary" spaces, but also in weighted spaces. For any m > 0, 
we define the Hilbert space L'^{m) by 

L\m) = {/gL2(R2)|||/||„<oo}, where 

1/2 



(^y^^(l + |eP)™|/(OPde) =Kf\2, 
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where b{^) = (1 + I^P)"^^^- If m > 1, then L'^{m) ^ L^(R^). In this case, we denote by 
Llijn) the closed subspace of L'^{m) given by 



Ll{m) = e L'^{m) j w{C)d^ = Q 



Weighted Sobolev spaces can be defined in a similar way. For instance, 

H^{m) = {we L^{m) \ dif G L'^{m) for i = 1,2} . 

As is well known, the operator C is the generator of a strongly continuous semigroup 
e'^^ in L'^{m), see Appendix ^ Note that, unlike e*^, the semigroup e"^^ does not commute 
with space derivatives. Indeed, since diC = (£ + |)c?j for i = 1,2, we have die'^^ = e^e'^^di 
for all r > 0. Using this remark and the fact that V ■ v = 0, we can rewrite equation ([T^) 
in integral form: 

w{t) = e^^wo - / e-3(--^)v ■ e^^-'^^{v{s)w{s)) ds , (19) 
Jo 

where wq = w{0) G L'^{m). The following lemma shows that the quadratic nonlinear term 
in ( ]T9| ) is bounded (hence smooth) in L^(m) if m > 0. 

Lemma 3.1 Fix m > and T > 0. Given wi,W2 G C'^([0, T], L^(m)), define 

R{t)= [ V ■e^^-'^^{vi{s)w2is))ds , 0<r<T, 
Jo 

where vi is obtained from Wi via the Biot-Savart law (It). Then R G C°([0, T], L^(m)), 
and there exists Cq = Co{m,T) > such that 

sup ||i?(r)||„ < Co( sup ||wi(r)||„)( sup ||w2(i")||m) • 

0<T<T ^0<T<T ' ^0<T<T ' 

Moreover, Co{m,T) — as T ^ 0. 

Proof: Choose q G (1,2) such that q > Using Proposition [A.5| (with N = p = 2), 

we obtain 

KR{r)\2<C r L_^|5"vi(.)«;2(s)|,ds, 

Jo a^T-syi 

where a(r) = 1 — e"^ and 6(^) = (1 + I^P)^''^. Moreover, using Holder's inequality and 
Lemma p.l| , we have 

I^^Vit^alg < |&"W2|2|Vl|^ < C||w2||mkl|9 < C||«;2||m|kl||m , 

2-9 

since L^(m) ^ L^{R?). Therefore, for all r G [0,T], we find 



/ ^^ds)f sup ||wi(s)|U) f sup ||w2(s)|U) 

Jo a{s)i ^ ^0<s<T ' ^0<s<T ^ 



mr)\u<c 

This concludes the proof of Lemma |0|. QED 
We now show that (|l^ has global solutions in L^(m) if m > 1. 
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Theorem 3.2 Suppose that Wq G L'^{m) for some m > 1. Then (I4) has a unique global 
solution w G C°([0, 00), L^(m)) with w{0) = Wq, and there exists Ci = Ci(||wo||m) > 
such that 

lh(T)|U<Ci, r>0. (20) 

Moreover, Ci(||wo||m) as ||wo||m 0. Finally, ifwo G Ll{m), then f^2U!{^,T) dC, = 
for all T >0, and lim ||w(r)||m = 0. 

Proof: Using Lemma ^TT] and a fixed point argument, it is easy to to show that, for any 
K > 0, there exists T = T{K) > such that, for all Wq G L'^{m) with ||ifo||m < K, 
equation (|19D has a unique local solution w G C°([0, T], L^(m)). This solution w{t) 
depends continuously on the initial data wq, uniformly in r G [0,T]. Moreover, T can be 
chosen so that ||w(r)||m < 2||wo||m for all r G [0,T]. Thus, to prove global existence, it is 
sufficient to show that any solution w G C°([0, T], L'^{m)) of (pj]) satisfies the bound (|20|) 
for some Ci > (independent of T.) 

Let Wq G L'^{m), T > 0, and assume that w G C°([0, T], L^(m)) is a solution of (p^ . 
Without loss of generality, we suppose that T >T = T(||wo||m)- Since L'^ijn) ^ Lp(R^) 
for all p G [1,2], there exists C > such that |w(r)|p < C||iL'(r)||m < 2C||wo||m for all 
TG [0,T]. By we alsohave \w{t)\p < Cpa{f)-^+^/P\wo\i < CHwolU for aU r G [T,T]. 
Thus, there exists C2 = C2(||wo||m) > such that |w(r)|p < C2 for all p G [1,2] and all 
r G [0,T]. Moreover, C2(||wo||m) — > as ||wo||m — ^ 0. To bound ||^|™'w(r)|2, we compute 



2dr 



R2 



= J l^l^'^iwAw + ^{$-V)w + w^ - w{^r ■ V)w} d^ . (21) 
Integrating by parts and using the fact that divv = 0, we can rewrite 

\^\^"'w{Aw)d^ = - J \^f"'\Vw\^d^ + 2m^ J {^{^""-Wd^ , 

J\^rw{^^-v)wd^ = iy'iep™(vv)^Me=^y"i^rv-(v^2)de 

= -m [ le'^-'i^ ■ ^r)w' d^ . 



Inserting these expressions into ([2l|), we find: 

+ 2m^J lep'^-Wde + m J iei"""'(^ ■v)w7Me . 

We next remark that, for all e > 0, there exists > such that 

|e|'™~Wde < e / |e|""«;'de + a /w^'de. 



J ■ v)w' d^\ < tj\i\^^w'di + C,\^^\'^jw^di. 



Note that, just as in of Lemma p?T| |v(r)|oo < C'l'W^lp I'W^lg " where l<p<2<g<oo 
and ^ + = We aheady know that |w(r)|p < C2 for all r G [0,T]. Choosing 



■'- ~ 8m' ? = 4 + and using (|T^ to bound |w(r)|q, we obtain 



where C3 = 6*3(771, ||tyo||m)- Thus, we see that for any 5 G (0, m — 1), there exists e > 
sufficiently small so that 



+ a(i + c3)(i + r-3/^^) I^Me. 



Integrating this inequality, we obtain that /|^p™'w^d,^ < C4 for all r G [0,T], where 
C4 = C4(m, ||wo||m) — >■ as ||wo||m — ^ 0. This proves that all solutions of (|T9|) in L^(m) 
are global and satisfy (^). Finally, if ^ -^o("^) ^^"^ "^o is the velocity field obtained 
from wq via the Biot-Savart law ([l6|) , then vq G L^(R^)^ by Corollary p.3| . In this case, 
we already observed that |w(r)|2 converges to zero as r — > +00, and so does ||^|™'w(r)|2 



by (p2p. This concludes the proof of Theorem 3.2. QED 



We now explain the motivation for introducing the scaling variables and discuss briefly 
how we will proceed. The spectrum of the operator C acting on L'^{m) consists of a 
sequence of eigenvalues ad = { — | | k = 0,1,2,...}, plus continuous spectrum ac = 
{A G C I 3fJA < — (^^Y^)}, see Appendix ^. Ignoring the continuous spectrum for a 
moment (since we can "push it out of the way" by choosing m appropriately) we choose 
coordinates for L'^{m) whose basis vectors are the eigenvectors of C These eigenvectors 
can be computed explicitly. Expressed in these coordinates (|l^ becomes an infinite set 
of ordinary differential equations with a very simple linear part and a quadratic nonlinear 
term. In the study of dynamical systems, powerful tools such as invariant manifold theory 
and the theory of normal forms have been developed to study such systems of equations. 
We will apply these tools to study ([T^). In particular, we will show that given any 
fi > 0, the long time behavior of solutions of (|l^) up to terms of order 0{e~^'^) (which 
corresponds to the behavior of solutions of (0) up to terms of order 0{t~^^^^^ ) is given by 



a finite system of ordinary differential equations which results from restricting (|lj) to a 
finite-dimensional, invariant manifold in L^(m). This manifold is tangent at the origin to 
the eigenspace spanned by the eigenvalues in ad with eigenvalues bigger than — /i. Using 
these manifolds we will see that we can, at least in principle, calculate the long-time 
asymptotics to any order. As an example we will show that in computing the asymptotics 
of the velocity field we obtain additional terms beyond those calculated in and [ pil 
due to the fact that we work with the vorticity formulation. We also exhibit an example 
that shows that in general the asymptotic behavior of the solution of (0) contains terms 
logarithmic in t. 

There are some technical difficulties which have to be overcome in order to implement 
this picture. The two most important ones are the presence of the continuous spectrum 
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Fig. 1. The spectrum of the hnear operator L in L'^{m), when m — A. 

of £, and the fact that the nonhnear term in (0) is not a smooth function from L^(m) 
to itself. We circumvent these problems by working not with the differential equation 
itself, but rather with the semiflow defined by it. For this semiflow we can easily develop 
estimates which allow us to apply the invariant manifold theorem of Chen, Hale and Tan 



In order to apply the invariant manifold theorem we consider the behavior of solutions 
of (|14D in a neighborhood of the fixed point w = 0. Since we wish for the moment to 
concentrate only on this neighborhood we cut-off the nonlinearity outside of a neighbor- 

> R be a bounded C°° function such 
l^^llm > 2. For any tq > 0, we define 



hood of size vq > 0. More precisely, let x : L^(m 



that x(^) = 1 if \\w\\m < 1 and x(^) = if 
Xroi'^) = x{'^/t"o) for all w G L'^{m). 



Remark 3.3 Such a cut-off function x ^o,y '^ot exist for all Banach spaces. However, 
for a Hilhert space like L'^{m) one can always find such a function. 



We now replace (|T3) by 



dT-w = Cw — (v ■ V)(Xr 



\W]W] 



(23) 



Note that (^31) and (|1J) are equivalent whenever ^ "'"o; but ( psf ) becomes a linear 

equation for > 2ro. Using Proposition |A.2| and the analogue of Lemma |3.1| , it is 



easy to show that all solutions of (^) in L^(m) are global if m > 0, and stay bounded if 
m > 1. In the sequel, we denote by $^0'™ the semiflow on L^(m) defined by (pB]), and by 
$™ the semiflow associated with (14). The following properties of will be useful: 

Proposition 3.4 Fix m > 1, ro > 0, and let <|)^0'''" he the semiflow on L'^{m) defined by 
1^^). If vq > is sufficiently small, then for t = 1 we can decompose 



ro,rn 



A + n 



(24) 



where A = is a bounded linear operator on LP'im), and TZ : L'^{m) — > L?'{m) is a C°° 
map satisfying TZ{0) = 0, DTZ{0) = 0. Moreover, TZ is globally Lipschitz, and there exists 
L > (independent of vq) such that Lip (7?.) < Lvq. 
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Proof: Let Wi,W2 G L'^{m), and define Wi(r) = ^^^°'"^Wi for i = 1,2 and < r < 1. Then 

Wiir) = e^'^Wi - re-5(--«)V ■ e(^-^)^(v,(s)x.„(«;.(s))«;,(s)) ds , t = 1,2 , 
Jo 

where Vj(r) is the velocity field obtained from Wj(r) via the Biot-Savart law (|TB|). Pro- 
ceeding as in the proof of Lemma we find that there exists C > 1 and K > such 
that 

sup \\Wi{T)-W2iT)\\m<C\\Wi-W2\\m + Kro SUp || Wl (t) - (^) || m • (25) 

0<T<1 0<T<1 

Thus, assuming Ktq < 1/2, we obtain H^^^'^ti?! — $^°''"iy2||m < 2C||wi — W2\\m for all 
T e [0, 1]. We now define 7^ = - A, so that 

= - /"'e-5a-)v.e(^-^)^(v,(r)x.oK(r))ti;i(r))dr , ^ = 1,2 . (26) 

In view of (|25|) , we have \\7l{wi)—7l{w2)\\m < -^'"o||'"^i~W2||m, where L = 2Cfr. Moreover, 
using (p6D, it is not difficult to show that TZ : L'^{m) L'^{m) is smooth and satisfies 
Dn{0) = 0. QED 

In the rest of this section, we fix A; G N and we assume that m > k + 2. As is 
shown in Appendix ^, the spectrum of A = e^ in L'^{m) is ct(A) = Sc U T^d, where 

= {e-"/2 I n = 0, 1, 2 . . .} and = {A G C I |A| < e-^"^)}. Since m - 1 > k, this 
means that A has at least k + 1 isolated eigenvalues Aq, Ai, . . . , A^, where \j = e~^^'^. Let 
Pk be the spectral projection onto the (finite-dimensional) subspace of L'^{m) spanned by 
the eigenvectors of A corresponding to the eigenvalues Aq, . . . , A^, and let Qk = 1 — Pk- 
Applying to the semiflow $^0'™ the invariant manifold theorem as stated in Chen, Hale, 
and Tan [Q, we obtain the following result: 

Theorem 3.5 Fix k&'N,m>k + 2, and choose 1^1,^2 G R such that | < /ii < 
< Let Pk,Qk be the spectral projections defined above, and let = PkL'^im), 

Es = QkL'^{m). Then, forrQ > sufficiently small, there exists a and globally Lipschitz 
map g : Ef. ^ Eg with g{0) = 0, Dg{0) = 0, such that the submanifold 

Wc = {wc + g{wc) I G E^] 

has the following properties: 

1. (Invariance) The restriction to of the semiflow •J)^"''" can be extended to a Lipschitz 
flow on Wc- In particular, $^0''"(lVc) = W^c for all t > 0, and for any Wq G Wc, there 
exists a unique negative semi-orbit {w{T)}r<o in Wc with w{0) = Wq. If {w{T)}r<o is a 
negative semi-orbit contained in Wc, then 

limsup -j— j- In ||t(;(r) ||m < /ii • (27) 

T— ►— 00 \T\ 

Conversely, if a negative semi-orbit of satisfies 

limsup j^ln ||ty(r)||m < /i2 , (28) 
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then it lies in W^- 

2. (Invariant Foliation) There is a continuous map h : L?'{m) x Eg ^ such that, for 
each w G Wc, h{w,Qk{w)) = Pk{w) and the manifold = {h{w,Ws) + Ws \ Wg E Es} 
passing through w satisfies $^0'™(M^) C M^^o.'"(^) and 




Moreover, h : L'^{m) x Eg ^ E^ is in the Eg direction. 

3. (Completeness) For every w G L'^{m), M^^nWc is exactly a single point. In particular, 
{M^}u;gVKc is a foliation of L'^ijn) overWc. 



Remark 3.6 Although the nonlinearity TZ is C°° , the invariant manifold Wc is in general 
not smooth. However, for any p > 1 such that kp < k + 1, one can prove that the map 



g : Ec ^ Eg is of class C if tq is sufficiently small. (See ji^/ , Section 6.1.) Note also 
that the manifold Wc is not unique, since it depends on the choice of the cut-off function 
X- This is the only source of non-uniqueness however - once the cutoff function has been 
fixed, the construction of ^ produces a unique global manifold. 

Proof: All we need to do is verify that hypotheses (H.1)-(H.4) of Theorem 1.1 in 
hold for the semiflow ^^O'*". Assumption (H.l) is satisfied because w <^^0'^w is globally 
Lipschitz, uniformly in r G [0, 1], see (0). Hypothesis (H.2) is nothing but the decom- 
position •l)^"''" = A + 7?. obtained in Proposition |3]^. Assumption (H.4) is a smallness 
condition on Lip(7?.), which is easily achieved by taking ro > sufficiently small. To verify 
(H.3), we recall that L^(m) = Ec® Eg, and we set Ac = PkAPk, Ag = QkAQk- Since 
ct(Ac) = {1, e~2 , . . . , e~2 1 and all eigenvalues of A^ are semisimple, it is clear that A,, has 
a bounded inverse, and that there exists Cc> I such that 

< Cc e-^^\\w\\m , for all j G N and all w E Ec . 

On the other hand, take e = 0ifm>fc + 2ore>0 arbitrarily small if m = /c + 2. By 
Proposition |A.2|, there exists Cg > 1 such that 



||A-^w||m < Cg e 2 , for all j G N and all w e Eg . 

These estimates are exactly what is required in (H.3). QED 



By Theorem |3^ , there exists ri > such that all solutions of (|T^) with ||w(0)||m < ri 



satisfy ||w(r)||m < tq for all r > 0. For such solutions the semiffow $™ defined by ([T 



coincides with •J)^"''", and thus the invariant manifolds in Theorem p.5| are also locally in- 
variant with respect to As a consequence, there exist a multitude of finite-dimensional 
invariant manifolds in the phase space of the vorticity equation (|T^), hence in the phase 
space of the Navier-Stokes equation. Furthermore, points 2 and 3 in Theorem give us 
an estimate of how rapidly solutions approach these invariant manifolds. Since this result 
will be crucial in the applications, we state it here as a Corollary. 
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Fig. 2. A picture illustrating the invariant manifold Wc and the foliation {M^}u,gvKc of L'^{m) over Wc- 
The leaf Mq of the foliation passing though the origin coincides with the strong-stable manifold Wl°^ in 
a neighborhood of the origin. Note that dim(i5c) < oo, whereas the subspace Eg is infinite-dimensional. 



Corollary 3.7 Fix k E N , m > k+2, and let Wc be the submanifold of L'^{m) constructed 
in Theorem W^. Define 



W^"" = Wcn{we L\m) I \\w\\^ < ro} . (30) 

Then Wjf"^ is locally invariant under the semiflow $^ defined by (||^. If {w{T)}r<o is a 
negative semi-orbit of such that ||w(r)||m < for all r < 0, then w{t) G Wl°'' for 



all T < 0. Moreover, for any < /i2 (where fj,2 is as in Theorem \3. 3i) , there exist r2 > 
and C > with the following property: for all wq G L'^{m) with ||wo||m < ^2, there exists 
a unique Wq G Wjf^ such that ^'^{wq) G Wjf'' for all t > and 

||<l>:"(*o) - ^riwo)\lm < Ce-^" , r > . (31) 



Proof: It follows immediately from point 1 in Theorem |3.5| that Wjf^ is locally invariant 



under the semiflow $™ and contains the negative semi-orbits that stay in a neighborhood 
of the origin. Choose ri > so that Hitiollm < '"i implies ||<l'™(wo)||m < '"o for all t > 0. By 
points 2 and 3 in Theorem |3]^, given wq G L'^{m), there exists a unique wq G M^q Pi Wc. 



Setting Wq = Wc + g{wc), we see from the definitions that Wc solves the equation 

Wc = h{wo,g{wc)) ■ (32) 

Since h is continuous and Wc ^— ^ h{wo, g{wc)) is Lipschitz (with a small Lipschitz con- 
stant), it is clear that ( p^ ) has a unique solution Wc which depends continuously on 
Wq. Moreover, by (|29|) , Wc = if Wq = 0. Therefore, by continuity, there exists 
r2 G (0,ri] such that, if ||wo||m < ^2, then ||wo||m = ll'W^c + fi'('W^c)||m < i^i- In this 
case, max(||<l>;"(wo)|U, ||$™(wo)|U) < tq for all r > 0, and (|TD follows from (||). QED 



In particular, given /i < /i2, it follows from Corollary |3]^ that a solution w(r) of (p^ ) 
on cannot converge to zero faster than e"^"^ as r — >• -|-oo, unless w{t) = 0. For this 
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reason, Wjf"^ is usually called the "weak-stable" manifold of the origin. In the applications, 
we will also be interested in solutions that approach the origin "rapidly", namely with a 



rate (9(e~^^'^) or faster. By Theorem p.5| , all such solutions lie in the leaf of the foliation 
{Mw}weWc which passes through the origin. We refer to this leaf Mq as the "strong- 
stable" manifold of the origin, and we denote by Wj"'^ its restriction to a neighborhood 
of zero. In contrast to Wjf"^, the local strong-stable manifold Wj."'^ is smooth and unique. 
In particular, it does not depend on the way in which the cut-off function x is chosen. 



Theorem 3.8 Fix kElSi,m>k + 2, and let Ec, Eg be as in Theorem \3. 4 Then there 
exists ra > and a unique C°° function f : {ws G Eg \ \\ws\\m < r^} — > Ec with /(O) = 0, 
Df{0) = 0, such that the submanifold 



W. 



loc 



{Ws + f{Ws) \ WsE Es 



satisfies, for any fi G 



'k fc+l l 
^2' 2 J 



W. 



loc 



w G L (m) 



\w\ 



< r3 , limsupr Mn < — [ 

r— >+oo ■' 



(33) 



where |||u^|||-m = max(||Pfcw||m, ||Qfcit'||m)- In particular, if Wq G W^"'^, there exists T > 
such that $™wo G 1^]°^ for allT>T. 

Proof: Choose > sufficiently small so that ||$™Wo||m < '"o for all r > whenever 
||wio||m < 2r3. Take the function h{-, ■) of point 2 in Theorem p.5|, and define 



f{ws) = h{0,Ws) , for all Wg G Eg with 



Then / is of class C^, /(O) = 0, and the characterization ( p3D with yU = /i2 follows 
immediately from (^9]) (with w = 0). In particular, / is unique. Moreover, since any 



solution w{t) on W^"'^ converges to zero as r 

M holds for any /i G ^ 



+00, and since ^2 G (|, ^) was arbitrary, 
2 7 5 ij-v^ij-v^v^ ^ — 2 also. 



it is clear that (p3D holds for any n G (|, ■^), hence for fi = also. Finally, the 
smoothness of / and the fact that -D/(0) = can be proved using the integral equation 
satisfied by /. (See Section 5.2.) QED 

If W^"'^ is as in Theorem R]^, we also define the global strong-stable manifold of the 
origin by 



{wo G 
{wo G 



<I)"(wo) G Wi"" for some r > 0} 
limsupr""'^ In ||w(r)||m, < — /i} , 



(34) 



where /i G (|,^]. Then W, is a smooth embedded submanifold of L'^{m) of finite codi- 



mension. (See Henry |17], Sections 6.1 and 7.3.) The proof of this statement is rather in- 



volved and uses two main ingredients. First, for any r > 0, the map : L 
is one-to-one. This property of the semifiow <I>™ is usually called backwards uniqueness. 
Next, since preserves the integral of w, it is clear that Ws is contained in the subspace 
{w G L\m) I J^, 



mm 



""^(0 = 0}. If w(r) is any solution of (|T4D in Ll{m) and if 



v(r) is the corresponding velocity field, it follows from Corollary [B.3| that v(r) G L^(R^) 
Then, the classical energy equality 



V r 



|v(0)| 



lVv( 



lids 



r > 
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shows that |v(r)|2 is a strictly decreasing function of r unless Vv(r) = 0, which is 
equivalent to w{t) = 0. In other words, the semiflow $™ is strictly gradient in the 
subspace Ll{m). 



4 The long-time asymptotics of solutions 

In this section we give three applications of the results of the preceding section. In the 
first, we examine the long-time asymptotics of small solutions of (JT^) and show that all 
such solutions with non-zero total vorticity asymptotically approach the Oseen vortex. 



thereby recovering the results of Giga and Kambe |Tj]. The invariant manifold approach 
yields additional information and as an example, we show that all solutions approach the 
vortex in a "universal" way. 

Prior investigations of the long-time asymptotics of the Navier-Stokes equations (P, 
||T0| [pO|] ) have yielded expressions in which the terms were proportional to inverse pow- 
ers of y/t. As a second application of the invariant manifold approach we show that if 
one extends these calculations to higher order (and in principle the invariant manifold 
approach allows one to extend the asymptotics to any order) one must include terms in 
the asymptotics proportional to (log(t))"/ (Vi)^, where a, /5 G N. We also exhibit specific 
classes of solutions for which such logarithmic terms appear. 

Finally, in a third application we extend some recent results of Miyakawa and Schonbek 
20| on solutions of the Navier-Stokes equations that decay "faster than expected" . The 
invariant manifold approach allows us both to give a more complete characterization of the 
set of such solutions and provides a natural geometrical interpretation of the conditions 



m 
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4.1 Stability of the Oseen Vortex 

We begin by considering the behavior of small solutions of (p!4D in the space L^(m) with 
m = 2. Acting on L'^{2), the operator C has a simple, isolated eigenvalue Aq = 0, with 
eigenfunction 

G(0 = 2^ e-l^l^/\ eeR^ 
An 

(See Appendix |B[) Let v*^ denote the corresponding velocity field, satisfying rotv"^ = G. 
From the explicit expression ( |101| ), it is clear that ^ _L v*^(,^), hence 

■ VG = , ^ G (35) 

As a consequence, for any a G R, the function w{^) = aG{^) is a stationary solution of 
(|l^), called Oseen's vortex. 

Let Ec = span{G}, let Eg be the spectral subspace of C corresponding to the contin- 
uous spectrum = {A G C | 3?(A) < — |}, and consider the local center manifold W^°'^ 
given by Corollary p.7| (with k = 0, m = 2). We claim that 

W^"" = {aG I a G R , |a|||G||2 < tq} . (36) 

Indeed, if |a|||G||2 < tq, then w{C,,t) = aG{^) is a solution of (|1^) such that ||w(r)||2 < tq 
for all r < 0, hence w{t) G by Corollary Since C {aG + g{aG) | a G R} 
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for some g : ^ E^, it follows that g = and that (^) holds. (Remark that, in this 
particular case, the local center manifold is unique.) Applying Corollary we conclude: 



Proposition 4.1 Fix < < |. There exist positive constants r2 and C such that, for 
any initial data wq with ||wo||2 < ^2, the solution w{-,t) of ( fl^j satisfies: 



\\w{;T)-AG{-)h<Ce-^^, r>0, (37) 

where A = wo(0 d^- 

Proof: If r2 > is sufficiently small, it follows from (|3l|) and ( ^6| ) that ( ^7|) holds 
for some A G R. Now, an important property of ( [T^ ) is the conservation of mass: if 

«(^) = /r2 ^(^' ^) then 



/ {Cw - V ■ Vw) d^= I V ■ (Vw + -^w - vw) d^ = . (3^ 

iR2 JRZ 2 



Since Jj^a G'(^) d^ = 1, the conservation law (|38|) implies that A = Jj^a wo{^) d^ in 
QED 



To facilitate the comparison of our results with those of we revert to the unsealed 
variables (x, t). Let 

Q^a:,t) = ^G(^^] , u^(x,t) = ^^v«^ ^ 



1 + t VyiTt/ VT+t \VTTt 

Thus is the solution of (0) corresponding, via the change of variables (0), to the solution 
G of (p!^ , and is the associated velocity field. From Proposition [4.1| , we obtain: 

Corollary 4.2 Fzx < yU < |. There exists r2 > sfic/i i/iat, /or a// initial data cjq ^ 
L^(2) with llcuolh < ^"2; solution ijj{x,t) of satisfies 

\Lu{-,t)-An{-,t)\p< , , l<p<2, t>0, (39) 

(l + t)^+^"p 

where A = J^2^o{x)dx. Ifu{x,t) is the velocity field obtained from Lj{x,t) via the Biot- 
Savart law (Qj, then 

|u(-,t)-Av^(-,t)|, < l<g<oo, t>0. (40) 

(1 + t)2+^-^ 

Proof: Let u;(x, t) be the solution of @) with Ci;(-, 0) = ujq, and let t(7(,^, r) be the solution 
of ([I|) with the same initial data. If 1 < p < 2, then L'^{2) ^ ^^(R^). Using (|T2|) and 
(^), we thus obtain 

\u;{;t)-AQ{;t)\, = {l + t)-'+-^\w{;log{l + t))~AG{-)\, 

< C(l + t)-^+-p\\wi-, log(l + t)) - AGi-)\\2 < Gil + ty'-^^'p . 
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It then follows from Lemma |2.1| that (^) holds for all q G (2, oo). Finally, assume 
that 1 < g < 2, and fix m G (2/g,2). If w{t) = w{t) — AG and if v(r) denotes the 
corresponding velocity field, it follows from Proposition [B.l| and Holder's inequality that 



|v(r)|, < C|6™-2v(r)|4 < C|6™«)(r)|2 < C||w)(r)||2 < Ce 



r > , 



where b{^) = (1 + Using the change of variables ([l3|), we thus obtain (^) for 

1 < g < 2. QED 



Remark 4.3 Once is known for p G [1,2], a bootstrap argument using the integral 
equation satisfied by uj{x,t) gives the same estimate for all p G [l,C)o] if t > 1 (see / [7^/ , 
Proposition 5.3.) Similarly, ^4^ ) holds for q G (l,+oo] ift > 1. However, as we shall see 
below, the difference u{t) — A\^{t) is not in L^(R^)^ in general. 



Remark 4.4 In Giga and Kambe show that ^3^) remains true if uq is any finite 

measure on R^ satisfying /j^2(l + |a^P)|t^o|(dx) < oo and |i^o|(dx) -C 1. Note that, in 
this case, the solution uj{x,t) satisfies \\uj{-,t)\\2 < r2 for all t > 1, so that Gorollary \4.3{ 
applies to uj{x,t) restricted to t > 1. In 1^, Garpio also shows that solutions of (|^ with 
finite measures as initial data satisfy estimate Instead of requiring ujq to be small 

in some norm, she assumes that the Gauchy problem for the vorticity equation with 
initial data A6 has a unique solution, where A = ^oi^x) and 6 is the Dirac mass at 
the origin. However, to the best of our knowledge, the uniqueness of solutions for this 
problem has only been established if \A\ is sufficiently small. 

We next note that not only does the invariant manifold approach yield the stability 
of the Oseen vortex but it also allows us to systematically compute the way in which the 
vortex is approached. Somewhat surprisingly, the approach to the vortex solution has a 
universal form if the norm of the initial data is not too large. More precisely, we prove: 

Theorem 4.5 Fix ^ < /i < 1. There exist r2 > and C > such that, for all initial 
data Wo G -^^^(3) with WwqWs < r2, the solution w{-,t) of (|7^ satisfies 

\\wi^,T) - AGiO + liB^^l + 526)G(0e-^/l3 < Ce''^^ , r > , (41) 
where A = J wo(0 d^; = / 6^^o(0 and B2 = J ^2Wo{^) d^. 

Proof: Acting on -^^^(3), the operator C has, in addition to the simple eigenvalue Aq = 0, a 
double eigenvalue Ai = — |, with eigenf unctions Fi{^) = — yG(^) and -^2(0 = 
(See Appendix 0.) Let Ec = span{G, Fi, F2}, let Eg be the spectral subspace of C 
corresponding to the continuous spectrum Uc = {A G C | 5R(A) < 1}, and consider the 
local invariant manifold Wj:"'^ given by Corollary (with k = 1, m = 3). 

By construction, a point w G Wj:"'^ can be written, in a unique way, as w = aG + 
PiFi + P2E2 + g{a,f3), where f3 = {^1,^2), and g{a,f3) G Eg. The coefficients a,f3 are 
given by the formulas 

«= / wiOd^, Pi = -[ 6^(0de, 2 = 1,2. (42) 
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Now, if w{-,t) evolves according to (0), we know from ( pH]) that a = 0, i.e. a does not 
change with time. In an analogous way, we find 



A = -/ U^w-vVw)d^ = -l/3i, 1 = 1,2. (43) 

Indeed, since V ■ v = and rot v = diV2 — <92f i = w, we have the identities 

^iCw + ^^iw = di {^idiw + ]^ilw -n?j+d2 {iid2W + , 
^^v ■ Vw = di {^iiviw - viv-^ + 82 {iiV2W + ^(^^1 - vl)^ , 
which prove (|^) for i = 1. The case i = 2 is similar. 

Remark 4.6 In the unsealed variables, ( j^^ means that equation ^ conserves the first 
moments of the vorticity uj{x,t). 

Thus, remarkably, the semiflow induced by ([llD on the three-dimensional invariant 
manifold is described by the linear equations (0), (0)' To complete the investigation 
of how the solutions evolve on the center manifold, we must estimate the nonlinear term 
g{a, f3). Given | < /i < 1, we claim that there exist e > and Cg > such that 

\\g{a,m3<CM' ^ (44) 

for all («,/3) € R'^ with |a| + \/3\ < e. Indeed, we know from Theorem (see also 
Remark ^ that ^ : R^ is of class C^m, and that c/(0,0) = 0, £'^(0,0) = 0. 

Moreover, since aG is a fixed point for (p^lf) , in must lie in Wjf^ when a is sufficiently 
small, so we have g{a,0) = for such values of a. Thus, to prove (0), it is sufficient to 
show that dj3g{a,0) = when |a| is small. 

To prove this, we linearize the map ^^"'"^ in (E^ about the fixed point a*G obtaining 



Da^G^?'"" = A + D^^gTI . (45) 

Since TZ is smooth, and TZ = 0, DqR. = 0, Da*GTZ can be made arbitrarily small in norm 
by taking a* sufficiently small. In particular, since we know the spectrum of A explicitly, 
we can choose a* sufficiently small so that the spectrum of D^.g^^"'™" is "close" to the 
spectrum of A. More precisely, if < 5 < min(i, 1 — /i), there exists ao > such that if 
< a* < tto; -Da'G^i"'™ has a three-dimensional spectral subspace E'^* with the spectrum 
of Da*G^^i'^\Ef consisting of eigenvalues with absolute value greater than exp(— (^ + 5)) 
(these are the perturbations of the eigenvalues 1 and e~^/^ of A) and the remainder of the 
spectrum of Z^q.g'^i"'™ contained in a disk in the complex plane of radius exp(— (1 — 5)). 
Applying the invariant manifold theorem of [0] to the semiflow $^0'"^ about the flxed 



point a*G (instead of the origin), we flnd just as in Theorem p.5| that <|)!|^0'™ has a three- 



dimensional invariant manifold , which in is tangent to E'^ at a*G. Furthermore, 



as we noted in Remark |3.6| , once the cutoff function Xtq in (Pg) is flxed, the invariant 
manifold constructed in |^ is unique. Therefore, it follows from the characterizations 
(p7|), (p8|) that W^* actually coincides with Wl°'^ in a neighborhood of a*G. We complete 
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the argument by showing that E^* = span{G, Fi, F2}. This imphes that di3g{a*,0) = 
for a* sufficiently small, thus proving 



To compute E" , note that if we linearize ([T^) about the fixed point w = a*G, we 
obtain 

drw = C''*w = A^w + -(I ■ V^)^ + w-a*i^v^ ■ V^w + v ■ V^G) , (46) 

where v*^ is the velocity field associated with the Oseen vortex G and v is the velocity field 
constructed from w via (0). Differentiating identity ( |35D with respect to (j = 1,2), 
we obtain 

= a,(v^ . VgG) = iidjY^) ■ V^)G + (v^ ■ V^)idjG) 

= (v^^- ■V^)G + Y^ ■ Vi:Fj . (47) 



Combining (|35D and (|4^) we see immediately that 

= CG-a*{^r^ ■Vi:G + ^r^ -V^G) =CG = 
L-'F, = CF,-a*{^r^-V^F,+^r^^.V^G) = CF, = -^F,, j = l,2. 

Exponentiating these results we find that Da*G^\°'"^ = exp(£'^*) has eigenvalues 1 and 
e~2 with eigenspaces {G} and {^1,^2} respectively. This conludes the proof of (|^). 

Assume now that w{-,t) is a solution of ( |14|) on Wj:"'', and let a(r) = A, Pi{t) = 
Bie~^'^, (32{t) = B2e-^/\ where A, Bi,B2 are as in Theorem W^. Then 



«;(-, r) = «(r)G + A(r)Fi + /32(r)F2 + (^(^(r), /3(r)) 



hence (^Tj) follows directly from (|4^). 



On the other hand, if ||wo||3 < T2 and if w{t) is the solution of ( p!^ with initial data wq, 
Corollary |3^ shows that there exists a solution w{t) on PV^^"'^ such that \\w{t) — w(r)||3 < 
(^g-MT for all r > 0. In view of the previous result, this means that (PT| ) holds for 
some A,Bi,B2 G R. But equations (|38|) and (^31), which hold for any solution of (p!^, 
imply that A = / 1^0(0 d.^ ^-^id Bi = J CiWo{^)d^, i = 1,2. This concludes the proof of 



Theorem 4.5. QED 



Remark 4.7 Note that although the preceding computation of the asymptotics of the so- 
lutions near the Oseen vortex applies only if a* is small, and — ^ are eigenvalues of C""* , 
with eigenfuctions G and Fi,F2, for all values of a* . 



One can also rewrite the result of Theorem [4.5| in terms of the unsealed variables as 
we did in Corollary [4.2| . Denote 

2 

A ^/ X \ \ ^ Bi ^ / X \ 



1=1 

A. { \ X ^ TP. / X 
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Corollary 4.8 Fix | < /i < 1. There exists r2 > such that, for all initial data loq G 
L^(3) with llwolls < T2, the solution uj{-,t) of ^1^ ) satisfies 

H-,t)-iu^ppi-,t)\p< , , 1<P<2, t>0, 

where A = J ljq^x) dx and Bi = J XiUJo{x) dx, i = 1,2. If u{x,t) is the velocity field 
obtained from uj{x,t) via the Biot-Savart law 1^, then 

C 

|u(-,t) - Uapp(-,t)|g < TT — 1 < g < oo , t>0. 

(l + t)2+^-^ 

4.2 Secular terms in the asymptotics 

We next show that, in contrast to prior investigations of the long-time asymptotics of 
solutions of (0) which yielded expansions in inverse powers of the time, one will in general 
encounter terms in the asymptotics with contain factors of log(t). These terms arise from 
resonances between the eigenvalues of the linear operator C, and as we will see, computing 
them is straightforward using ideas from the theory of dynamical systems. 

In the previous subsection we saw that, if wq G L^{2) is sufficiently small, the solution 
w{C,,t) of (|14D with initial data wq approaches the Oseen vortex AG{^) as time goes to 
infinity, where A = J Wo(0 ^^'^ examine in more detail what happens ii A = 0, 

namely Wq G Ll{m). In this invariant subspace, we know from Theorem that all 
solutions converge to zero as r — +oo, without any restriction on the size of the initial 
data. Assuming that Wq G i^^(4), we will compute the long-time asymptotics of the 
solution up to terms of order 0{e~'^'^), for any fi in the interval (1, |). 

We first note that the operator C acting on -Lq(4) has, in addition to Ai = — |, a triple 
eigenvalue Ai = —1, with eigenf unctions 

m) = \m"- 4)G(o , H,{o = \{ei - e2)G{o , ^^3(0 = ^^^^(o . (48) 

Of course, these eigenf unctions are not unique, but the choices above are convenient ones 
for computation. 

Let Ec = spa.n{Fi, F2, Hi, H2, H3}, and let Es C Ivq(4) be the spectral subspace of C 
corresponding to the continuous spectrum cTc = {A G C | 3?(A) < — §}• Any function 
w G Lq{4) can be written as 

2 3 

wiO = 5^ Ai^.(0 + J2^jH,{0 + RiO , (49) 
i=i j=i 

where R E Eg. The velocity field v(,^) associated to w has a similar decomposition: 

2 3 

v(0 = E/^^^'''(0 + E^^-^''^(0 + v^(0 , (50) 

i=l j=l 

see Appendix 0. The coefficients Pi are given by (|4^) , and the corresponding formulas 
for 7j read: 

7,= / PMMOd^, J = 1,2,3, (51) 



21 



where p,{0 = - 4), p,iO = - ^i) and p3(0 = ^16- 

Assume now that w(^,r) is a solution of (0) in and consider the evolution 

equations for the coefficients 7j and the remainder R in (^91). Proceeding as in P3|), we 
find 

ii = -l3-l Pi(0(v- V)w7de, J = 1,2,3. (52) 
The following elementary result will be useful: 

Lemma 4.9 Assume thatw G Lq(2) nif^(2), and let\ be the velocity field obtained from 
w via the Biot-Savart law ^Tdj). Then, for any quadratic polynomial p{^) , 

[ p(0(v ■ v)w de = / {viV2{dip - dip) - {vi - vi)d,d2P) de . (53) 

JR2 JR2 

Proof: Since V ■ v = and w = diV2 — d2Vi, we have the identity 

p(v ■ V)w = viV2{dlp - dip) - {vl - vi)did2P + diEi + ^2-^2 , 

where Ei = pviw — ViV2dip + \{v'l — vi)d2P and E2 = pv2W + viV2d2P + \{y1 — vi)dip. 
Applying Proposition |B.1| with m = 2, we see from (|110|) that (1 + |^P)v G L°^(R^)^, 
hence E^ G L^(R^) for i = 1,2. In addition, since Vv G L^(R^)^ by Lemma |2.1| , we have 
diEi G L1(R2) for i = 1, 2, and (0) follows. QED 

It follows in particular from Lemma |4.9| that PiiOi^ ' V)wd^ = 0. Thus, surpris- 
ingly enough, the equation for 71 is linear: 71 = —71. In contrast, the equations for 72, 
73, and R contain nonlinear terms. For the purposes of this subsection, it will be useful 
to write out separately the terms that are quadratic in /3 = (/5i,/?2)- A direct calculation 
shows that 

(/3iv^^ + /52v^^) ■ v(/3iFi + P2F2) = m - f3i)^i^2 - Pimi - Q)m , (54) 

where $(0 = (87r2|^|4)-ie-l«l'/4(e-l«lV4 _ i + \^\y4). Thus, defining 

" ^ L ^'^'^^^^ ^ ^ Ij^' " ^'^'^^^^ = 3^ ' ^^^^ 

we see that the quadratic terms in the equations for 72, 73 are respectively HP1P2 and 
—n{l3f — /3|). Therefore, the equations for f3,j,R (where 7 = (71,72,72)) have the 
following form: 

Pi = ~Pi , $2 = ~P2 , 71 = -71 , 

72 = -72 + Kf3,[32 + /2(/3, 7, R) , (56) 

73 = -73-<Pi-Pi) + f3{f3,J,R), 

Rr = £i? + /3i/32^2-(/3i -/52')^3 + V-F(/3,7,i?)-v^- Vi?, 

where \E'2 = (^1 ~ ^2)^ ~ and \E'3 = .^1^2'^ ^ '^H'^. Moreover, since we have written 
out explicitly the quadratic terms in /3, the remainder terms F in (|56D satisfy the 
estimates 

l/.l + ||F||4 < cm + |7|)(|7l + Ili^lU) + C\\R\\l . (57) 
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To compute the long-time asymptotics of the system (p6D, a natural idea is first to 
study the behavior of the solutions on the five-dimensional invariant manifold Wj:"'^ tangent 
to Ec at the origin, and then to use the fact that all solutions of ( ^6] ) approach Wjf"^ faster 
than (9(e~'^^) for any /i < |, see Corollary |3.7| . This is the approach we follow, but we 
begin with a simplifying observation. Note that if one drops the terms V ■ F — ■ V-R 
from the last equation in (|56|) , the terms /3i/32^2 — {Pi — /^D^s ^-ct as a simple "forcing 
term" in the equation for R. As a result, one can find an explicit invariant manifold for 
the simplified equations, namely: 

R = g{(3) = -i5i(5^{L + \)-^^^ + - + \)-^^^ . (58) 

(Remark that \1'2, \I^3 G E^ and that the restriction of £ + 1 to is invertible.) As we will 
see, this manifold is a good approximation to an invariant manifold for the full equations 
(^). We thus introduce the new variable 

p = R- g{l3) , (59) 

and find that (|5B[) can be rewritten as 



1 • 1 

I3i = --A , 02 = --(32 , 7i = -7i , 

72 = -72 + «:/5i/52 + /2(/3, 7, P) , (60) 

73 = -i,-K{(3l-(3l) + UPn.p) . 
Pr = /:p + V-F(/3,7,p)-v^-Vp, 

where fj and F obey the same estimates as /-,■ and F, but with ||p||4 -|- |/3p substituted 
for 1 1 -R 1 1 4- Note further that, if R & Eg, then p E Eg as well. 



If we now apply the results of to the system (0), we find just as in Theorem ^ 
that (^) has a five-dimensional invariant manifold tangent at the origin to Eg. The next 
proposition provides an estimate of the behavior of this manifold. Let X be the Banach 
space X X Eg, equipped with the norm ||(/3, 7, i?)||x = |/3| + I7I + ||-R||4, and set 

= R^ X R^, with coordinates (/3, 7). 

Proposition 4.10 Fix p G (1, |) and 5 G (0, | — p.). There exists a map Q : Ec ^ Eg 
satisfying ^(0) = 0, DQ{Qi) = 0, such that in a neighborhood of the origin the graph of 
Q is left invariant by the semifiow defined by ij{6Q). Moreover, there exists r4 > and 
C > such that, for any solution {(3{t),'^{t), pij)) of (j^^ with initial data (/3g,7o,Po) 
satisfying ||(/3o5 7o'Po)||x < ^a, there exists a solution (/3'^(r), 7'^(r), ^(/3'^(r), 7'^(r)) on the 
invariant manifold with 

|/3(r) -/3^(r)| + \l{r)-Y{r)\ + l|p(r) - ^(^(r), 7^(r))||4 < Ce"'^^ . (61) 
Finally, there exists C5 > such that 

||^(/3,7)l|4<C,(|/3|='-^ + |7lH ■ (62) 
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With the exception of the last estimate, the proof of this statement is almost iden- 
tical to the proof of Theorem |3]^ and Corollary We prove (|62D in Appendix and 
concentrate here on showing how it can be used to determine the asymptotics of (|60D. 

We begin by studying the behavior of the solutions on the invariant manifold. Such 
solutions satisfy the system of ordinary differential equations: 

$1 = , $2 = ~(^2 , 7i = -7i , 

72 = -72 + «:/?i/32 + /2(/3,7,6^(A7)) (63) 

73 = -iz-^{Pl-Pl) + U(i,iM(in)) ■ 

There are various ways to analyze the asymptotics of solutions of (|63|) , but perhaps the 
simplest is to make the change of variables: 

Ti = 71 , = 72 + ^^Plf^2 \ log |/3i/32 1 1 , (64) 
Ts = 73 - n{Pl\\og{Pl)\ - \og{(il)\) . 

From a dynamical systems point of view, this is just a normal form transformation that 
eliminates the resonant terms in ([63|). Rewriting (|6^) in term of /3 and F = (ri,r2,r3), 
we find 

A = -^A , /32 = Af2 , Ti = -Fi , 

r2 = -r2 + ^2(/3,r) , Tg = -r3 + ^3(/3,r) , (65) 

where ^j(/3, F) is just fj{(3,j,Q{f3,^)) rewritten in terms of (/3, F) instead of {f3,^). 

Remark 4.11 The change of variables {f3,^) —>■ (/3, F) is not Lipshitz, so it might seem 
as if the standard existence and uniqueness theory for solutions does not apply to ^6di). 
However, due to the very simple form of the first equations in (65), we can first solve the 
equations for Pi and P2 explicitly and then insert these expressions into the equations for 
F. The equations for F are then non- autonomous, but Lipshitz in T and hence standard 
theorems imply that solutions of j^Udj ) exist and are unique. 

It is clear from (|65|) that /3i{T) = fojC"^''^ for i = 1,2 and Ti{t) = cie~^, where 
bi = A(0) cind Ci = ri(0). To determine the long-time behavior of r2 and T^, we note 
that e^Tjir) = Tj{0) + J^^ e'gj{f3{s),T{s)) ds. Thus, defining 

POD 

Cj = r,(0) + / e^^,(/3(s), F(s)) ds, j = 2, 3 , (66) 

and using the estimates on gj which come from (0), we find that for any /i G (1,3/2) 
there exists C > such that |rj(r) — e~'^Cj\ < Ce~'^^, for j = 2,3. Inverting the change 
of coordinates (Q), we immediately find: 

Lemma 4.12 Fix /i G (1, 3/2). There exist > and C > such that, for any solution 
of (^^ ) with initial data satisfying |/3(0)| + |7(0)| < r^, there exist constants 61, 62, c'^, 
and C3 such that, for all r > 0, 

A(r) = 6ie--/2 , P,iT) = 626-^/2 ^ ^^^-r ^ 

|72(r) -C2e-" + /s:re-^6i62|log|(6i62)|| < Ce"^" , (67) 
|73(r) - C3e-- - nre-^ibll log(6?)| - bj] log{bl)\)\ < Ce'^^ . 
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Remark 4.13 The constants c'j (j = 1,2,3) in (62) are related to Cj by 



Cl 



4 = C2 - log 1 6162 1 1 , 4 = C3 + K{bl\ log(6i)| - 62 1 log(^2)l) 



We now return to the full system of equations (|60D. We know that any solution 
(/3(r), 7(r), p(r)) of (|60|) with initial data in a sufficiently small neighborhood of the 
origin approaches a solution (/3'^(r), 7'^(r), ^(/3'^(r), 7'^(r)) on the invariant manifold with 
a rate 0{e~'^'^), see (0). On the other hand, from Lemma [4.12| , we know that there exist 
constants bi,c'j such that estimates (|67|) hold with fB'',^'^ substituted for (3,j. Finally, 
combining (|67| ) with (|62D, we obtain the bound ||^(/3'^(r), 7'^(r))||4 < Ce^^'^. Thus, the 



(68) 



following estimate holds for all solutions of ( |60D in a neighborhood of the origin: 

|/3i(r) - 6ie-^/2| + ^r) - 626-^/^! + \j,{t) - 4e-^| 
+ |72(r) -4e-" + Kre~"6i62|log|(&i62)|| + 



+ |73(r) - 4e-- - KTe~^{bi\ log(6f)| - log(6; 
Finally, we undo the change of variables (|59|). The result is: 



< Ce 



Theorem 4.14 Fix 1 < yU < |. If w{^,t) is any solution of ( |7^ in L'^{'^), then there 



exist constants bi, dp and C such that \\w{-,t) 



■Wappl-,i-j||4 



< Ce ^"^ for all r > 0, where 



app 



(e,r) 



e-n6iFi(0 + 62i^2(0) + + 4^^2(0 + 4^3(0) 

+ KTe~^{-bib2H2{0 + {bl - bl)Hs{0) (69) 

+ e--{-bM^ + l)-'^2(e) + {bl - bl){C + l)-^^3(0) • 



Proof: Let w{i,T) be a solution of ([I3) in Ll{A), and define /3(r), 7(r), /?(r) by (|Sp and 
p(r) by (^%. Since ||w(^, r)||4 — as r — +00, we can assume without loss of generality 
that /3(0), 7(0), p(0) satisfy the assumptions of Proposition [4.10| and Lemma |4.12| . In 
particular, (|68|) holds for some constants 6j, G R. The estimate ||w(-, r) 



appV 1 



< 



Ce ^'^ is then a direct consequence of (^), (|58|), and (|68D . QED 

The terms in (^) we wish to call particular attention to are those proportional to 
re^"^. When we revert to the original, unsealed, variables {x,t), these yield terms of the 
form (l + t)~^ log(l + t), which should be contrasted with previous asymptotic expansions 
which yielded only inverse powers of t. Proceeding as in the proof of Corollary ^4.2| , we 
obtain the following result in the original variables: 

Corollary 4.15 Under the assumptions of Theorem \4 .141 , the solution uj{x,t) of with 
initial data a;(x, 0) = wo{x) satisfies 



\uj{-,t) - UJ^pp{-,t)\p < 



Cn 



+ 



1 < p < 2 , t > 



[x,t) = j^^W, 
d from 

|u(-,t) - u 



Tog(l+t)). Similarly, z/ u(x, t), Uapp(x, t) are the velocity 



where LOgjppya,, V ) — j^u^appv ^i+t ' 

fields obtained from uj{x,t),ujs,pp{x,t) via the Biot-Savart law then 



appV ) 



l<g<cx), t>0. 
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Remark 4.16 Two other references which examine the asymptotics of solutions of (j^ to 
this order are and / (7^/ . We provide a detailed comparison of our asymptotics with these 
two references in /|7^/, but just mention here that since both of these references require that 
the initial velocity field satisfy (1 + |a;|)uo G L^(R^)^ in order to derive the asymptotics, 
the coefficients hi and 62 o^^e both zero for all the solutions they study. (This follows 
from Corollary \B.4 .) Thus many solutions of finite energy (i.e. of finite norm) are 
excluded from consideration by this hypothesis, and in particular one does not observe the 
logarithmic terms in the asymptotics. This is a further reason that we feel it is preferable 
to impose the decay conditions on the vorticity field rather than on the velocity. 



4.3 Optimal Decay Rates 

It has been known for a long time that there is a relationship between the spatial and 
temporal decay rates of solutions of the Navier-Stokes equation in R^, N > 2. For 
instance, for any initial data uq G L^(R^)^ fl L-'^(R^)^, there exists a global weak 
solution of the Navier-Stokes equation satisfying 

|u(t)|2<C(l + t)-^/\ t>0, (70) 

2l, p|. If, in addition, (1 + |x|)uo G L\R^)^, if follows from Wiegner's result [||] 



see 
that 



|u(t)|2 < C(l + t)-(^+2)/4 ^ ^ > Q _ ^7;^^ 

In [20|, T. Miyakawa and M.E. Schonbek investigate the optimality of the decay rate (|7lD . 



More specifically, they prove: 



Theorem 4.17 ^ Assume thatuo G ^^(R^)^, V-uq = 0, and (l + |a;|)uo G ^^(R^)^, 
N >2. Let u{t) be a global weak solution of the Navier-Stokes equation with initial data 
Uq satisfying (f7j|/ For all k,i E {1, . . . , A^}, define 



5m 



Then 



POO P 

Xi{uo)kix)dx , Cke= / Uk{x,t)ui{x,t)dxdt . (72) 
Rjv Jo Jn^ 



lim t^|u(t)|2 = (73) 



t— >oo 

if and only if there exists c > such that 



bki = and Cki = c5m , kJe{l,...,N}. (74) 



The proofs in are clear, but do not provide much intuition as to the meaning of 
the conditions (|7^ . As the authors themselves remark, "We know nothing about the 
characterization of solutions satisfying (cfc^) = (c(5fc^)." As we will demonstrate below 
in the case N = 2 (see also [|13| for the three-dimensional case), the invariant manifold 



approach provides a simple geometrical explanation of the meaning of both conditions in 
(^). In addition, it will allow us to construct additional solutions satisfying (|7^) , but 
which do not fit into the setting of Theorem |4.17| , because (1 + |x|)uo ^ L^(R^)^. Note 



however that our work does require that we consider solutions of ([T^) in the space i^^(4). 
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and this decay requirement on the vorticity does not appear in the work of Miyakawa and 
Schonbek. 

As before, we work with the vorticity equation rather than the Navier-Stokes equations 
themselves. This seems to us a particularly natural choice when examining the long- 
time behavior of solutions since, as we have noted above, the long-time asymptotics 
are influenced by assumptions about the spatial decay of the initial conditions, and such 
hypotheses are preserved by the evolution of the vorticity equation (see Theorem |3.2D , but 
not by the Navier-Stokes evolution. Since our goal is to recover the results by Miyakawa 
and Schonbek, we restrict ourselves to velocity fields that belong to L^(R^). At the level 
of the vorticity, this condition is equivalent to 



w{Od^ = 0, / ^MOd^ = 0, / 6^(0 de = 0, (75) 

see Corollary p.4 Thus, we shall study the solutions of the vorticity equation ([T^) in the 
invariant subspace of i^^(4) defined by (|75D. If w{^,t) is such a solution and if v(^,r) is 
the velocity field obtained from w{^, r) via the Biot-Savart law (^, then w and v can be 
decomposed as follows 

w{^,t) = 7i(r)i/i(0+72(r)if2(0+73(r)i/3(0 + ^(e,r) , (76) 
v(e,r) = 7i(r)v^i(0+72(r)v^H0+73(r)v^^(0+v^(e,r) , 

where the coefficients •yj (j = 1, 2, 3) are defined by (|5T|). Setting /3 = in (|56|), we obtain 
the evolution system 

71 = -7i , 

72 = -72 + /2(0,7,^) , (77) 

73 = -73 + /3(0,7,^) , 

Rr = /:/? + V-F(0,7,i?)-v^-Vi?. 



Also, using ( ^21) and Lemma [4.9| , we see that 



/2(0,7,i?) = - / v,V2d^, /3(0,7,i?)=/ {vf-v'^)d^^ (78) 

JR2 JR2 

where v is given by (|76D . 

Now, let Ec = span{Hi, H2, H3}, and let Eg C i^^(4) be the spectral subspace of C 
corresponding to the continuous spectrum ctc = {A G C | 3?(A) < — §}• Given 1 < < |, 
Theorem |3.8| (with k = 2 and m = 4) shows that, for sufficiently small ri > 0, the set 
Wg""^ defined by (^3]) is a smooth (infinite-dimensional) manifold which is tangent to Eg 
at the origin. Since the equations for a and Pi are linear, it is clear that W^"'^ is contained 
in the subspace of i^^(4) defined by (^), and so is the global strong-stable manifold Wg 
defined by (|3^ . The following result gives various characterizations of Wg . 



Proposition 4.18 Fix 1 < /i < |, and assume that wo G L (A) satisfies ^73^). Let w{^, r) 
be the solution of ^1^ ) with initial data wq, and let v(^,r) be the velocity field obtained 
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from w{C,, t) via the Biot-Savart law (Qj. Decompose w{C,, t) according to (7t), and define 
the coefficients cm by (fT^j, where 



u(a;, t) 



X 



,log(l+t) . 



(79) 



Then the following statements are equivalent: 

1) Wq lies in the strong stable manifold Wg. 

2) lim e^\\w{-,r)\\i = 0. 

r— >oo 

3) lim t\u{-,t)\2 = 0. 

4) 7i(0) = 0, 72(0) = C12, and 73(0) = C22 - Cn. 



Remark 4.19 Note that the global nature of the strong stable manifold means that Propo- 



sition U-i^ is not limited to solutions of small norm, but applies to solutions of arbitrary 
size. 



Proof: Applying Theorem f4.14| with 61 = 62 = 0, we obtain 



^0) 



for some Ci, 02,03 G R. We shall show that statements 1), 2), 3), 4) in Proposition [4.18 
are all equivalent to 
5) ci = C2 = C3 = 0. 



Indeed, since the functions Hi, H2, H3 are linearly independent, it is clear from (|33|) and 
(13) that 1) ^ 2) ^ 5). On the other hand, it follows from (|7|) and (^ that the 
coefficients cj in (^OD are given by the formulas 

POO 

ci = 7i(0), c,=7,(0)+/ eV,(0,7(r),^(r))dr , j = 2, 3 , 

Jo 

see also (|66|). Furthermore, using ([78|) and the change of variables (^), it is straightfor- 
ward to verify that 

roc p 

C2 = 72(0)-/ eM 1^1 (e, r)t;2(e, r) de dr = 72(0)- C12 , 

Jo ^R2 

C3 = 73(0)+/ eM (t;i(e,r)2-t;2(e,r)2)dedr = 72(0) + cn-C22 , 



where Cki is defined in (^). Therefore, 4) <^=> 5). Finally, from Corollary [4.15| , we have 
|u(-,t) -Uapp(-,t)|2 < C(l + t)-^, where 



1 ^ 

^ ' 7=1 



X 



Clearly, |uapp(-,t)|2 = j^, where = if and only if ci = C2 = C3 = 0. Thus, 3) -v^ 5). 
QED 
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We are now able to give an alternative proof of Theorem 4.17| in the particular case 



where N = 2 and Wq = (9i(uo)2 — ^^2(^0)1 € -^'^(4)- Indeed, let w{^,t) be the solution 
of ( p^ with initial data Wq. Since uq G L^(R^)^, it follows from Corollary [B.4| that ([75| ) 
holds, hence w{^,t) can be decomposed according to (|76D . Moreover, by Corollary p.5| , 
the assumption (1 + |x|)uo G L^(R^)^ is equivalent to 72(0) = 73(0) = 0. 

a) Assume first that (^) holds, namely Wq G Wg. Then point 4) in Proposition [4.18 



shows that C12 = and cn = C22, hence the matrix (cke) is scalar. In addition, 71 (0) = 0, 
hence bke = by Corollary [B.5| . 

b) Conversely, assume that (|7^ holds. Then 71 (0) = by Corollary p.5| , and since 



72(0) = = C12, 73(0) = = C22 — Cii it follows from Proposition [4.18| that Wq G Wg. This 
concludes the proof. QED 

Note that as a consequence of these investigations, we see that there are many other 
solutions u(a;, t) of the Navier-Stokes equation that lie in the strong stable manifold (and 
hence satisfy the decay estimate (|73|) ), but which do not satisfy the moment condition 
(H-|x|)uo G L^(R^)^. Put another way, the conditions ( |7^ are the necessary and sufficient 
conditions for (ff^) to hold within the class of solutions satisfying (1 + |a;|)uo G L^(R^)^, 
but if one looks at larger classes of initial conditions then one finds that (^4|) is no longer 
necessary, and that a better characterization of solutions satisfying (|73D is that they lie in 
the strong stable manifold Wg. Indeed, as we observed above, the velocity fields and 
v^^ corresponding to the vorticities H2 and do not satisfy (1 + G L^(R^)^. So 

if we choose as an initial condition any point in the strong stable manifold Wg which does 
not lie in the hyperplane Eg of functions of the form (^) with 71 (0) = 72(0) = 73(0) = 0, 
we obtain a solution of the Navier-Stokes equation satisfying (^), but not ([74|). 

Of course, this argument requires that there be some point in Wg for which either 
72 or 73 is nonzero. (Proposition [4.18| implies that all points in must have 71 = 
0.) This could fail to occur only if Wg coincided with the hyperplane Eg, namely if Eg 
was invariant under the semifiow defined by (0). This in turn would happen only if 
/2(0,0,i?) = /3(0,0,i?) = for all R G Eg. We now demonstrate explicitly that this is 
not the case. Choose 

K(0=9ii7i(0=ei(l-|eiV8)G(0. 
Then CK = —\K, so that K E Eg. The velocity field corresponding to K is 



v«({) = ftv«.K) = ^U._2 



and a direct calculation gives: 



/3(0, 0,K)= I {{v^y - (v^f) de = ^ . 



R2 



647r 



This shows that Wg (f_ Eg. On the other hand, it is easy to see that any w E Eg which 
is radially symmetric actually lies in Wg, because (0) reduces to the heat equation for 
radially symmetric vorticities. Thus the intersection Wg fl Eg contains at least an infinite- 
dimensional subspace of L'^{m). 
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Fig. 3. A schematic picture of the dynamics defined by equations (|77| ) in a neighborhood of the 
origin. All trajectories approach a three-dimensional invariant manifold Wc at a rate C'(e^'''^) or faster. 
The strong-stable manifold Ws contains all solutions that converge to the origin faster than e~'^. The 
intersection Ws H Es , which is depicted here by two points only, contains in fact an infinite-dimensional 
subspace of L^{m). 



Summing up, we see that for solutions of the Navier-Stokes equation in R^, whose 
vorticity hes in i^^(4), the stong-stable manifold identifies exactly those which converge 
to zero faster than C/t. If the velocity field at t = satisfies (1 -|- |x|)u G L^(R?y, this 
gives a natural geometrical interpretation of the moment conditions of It also shows, 
however, that there are additional solutions which decay with a rate faster than C/t but 
which do not satisfy the decay condition on the initial velocity field. 

We conclude this section by describing a somewhat surprising property of the strong- 
stable manifold. Although it was constructed to be invariant with respect to the dynamics 
of the vorticity equation in the rescaled variables, it is also invariant with respect to the 
dynamics expressed in the unsealed variables. Let "^"^ be the semifiow on L'^{m) defined 

by m- 

Proposition 4.20 The semifiow ^^J" leaves the manifold Wg invariant. 



Remark 4.21 This should be contrasted with the behavior of Wj."'^. In order to obtain 
a manifold corresponding to Wj."'^ that is invariant with respect to it is necessary to 
go the extended phase space L'^{m) x R_,_, where R+ represents the time axis. (See the 
discussion in Section 3 of in particular Theorem 3.6, where this question is discussed 
in a related context.) 

Proof: The proposition follows from the characterization of Ws in point 3 of Proposi- 



tion^]!^. Since u(-, t) is the velocity field corresponding to the vorticity field '^'^{wq), and 
since this condition is clearly invariant with respect to time translation, the invariance 
of Ws with respect to \E']" is immediate. The fact that Ws is invariant with respect to 
both and is related to the fact that solutions in Wg can be characterized by their 



30 



asymptotics as t ^ +00. To be more explicit, note that combining point 2 in Proposi- 
tion |4.18| , with point 5 in the proof of that proposition, we see that Wq lies in Wg if and 



only if the solution w{-,t) of (JT^) with initial data Wq satisfies 

lim e^|w(-,r)|2 = . (81) 

T— >+00 

But letting u!{x,t) = j^w{^^;=,log{l + 1)), we see that (PI) is equivalent to 

lim t^/^\uj{-,t)\2 = . (82) 

t—r + OO 

Note also that lj{x, 0) = wq{x). Thus, Ws can also be characterized as the set of points 
in L^(4) satisfying ( fTSD and for which 

lim t3/2|^-(^o)|2 = . 

t— >+oo 

Since this characterization is again invariant with respect to \E'™(wo), this shows that Wg 
is an invariant manifold in the original variables as well as in the rescaled variables. QED 

A Spectrum of the linear operator 

In this section, we assume that G N, > 1. We consider the linear operator C given 
by 

C = A^ + y-V^ + ^ , ^eR"". (83) 
As in Section ^, we shall work in the weighted space L'^{m) defined by 

L\m) = {/gL2(R^)|||/|U<oo}, (84) 

Since its coefficients depend linearly on the space variable ^, the operator C becomes a 
first order differential operator when expressed in the Fourier variable p. Our convention 
for Fourier transformation is 

fip) = I /(Oexp(-ip-Ode, (85) 

/(O = ttA^ / fiP) exp(ip ■ ^) dp . (86) 

The expression of the operator C in Fourier space is: 

(A/)(p) = -(br + ^P-v,)/(p). (87) 

The aim of this section is to prove the following result, which underlies our approach 
for computing the long-time asymptotics of the vorticity equation: 
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Theorem A.l Fix m >0, and let C be the linear operator (85) in L'^{m), defined on its 
maximal domain. Then the spectrum of C is 

a{C) = {\ec\^{\)<^-'^]u{-\\ke^] . 

Moreover, if m > ^ ^.''^d if k G N satisfies k + ^ < m, then = —\ is an isolated 
eigenvalue of C, with multiplicity (^"^^ • 



In the one-dimensional case, this result is proved in |T^, Appendix A. We give here 
a slightly different proof, which is valid for all A > 1. We begin with a few elementary 
observations. 

(1) The discrete spectrum of C Fix A; G N, and take a = (ai, . . . , on) G such that 
|a| = ai + . . . + aN = k. Then the Hermite function (pa £ iS(R'^) defined by 

= (ip)"e-|Pl' =il"lK\..p7e-l^'l' , or 

1 

(47r)^/2 

is an eigenfunction of C with eigenvalue — |. Thus, for any m > 0, we have c^C) D 
{ — 1 1 k E N}, and the multiplicity of the eigenvalue = — | is greater or equal to 
rT^)=#{«GN^||a| = fc}. 

(2) The "continuous" spectrum of C. Fix A G C such that 3?(A) < A/4 and —A ^ N. 
The function ipx : — > R defined (in Fourier variables) by 



V^A(p) = b|-'"e-lPl , PGR 



is then an eigenfunction of C with eigenvalue A. It is clear that ipx G C°°(R^), and a 



standard calculation (see for instance |jT6[, Section 2.3.3) shows that 



22^7rfr(A) 

In particular, G L'^{m) if and only if 3?(A) < ^ ~ T- Since the spectrum of C is closed, 
this shows that (t(A) D {A G C 1 3fJ(A) < f - f }. 

(3) The spectral projections. For all a = (ai, . . . , a^v) G N^, we define the Hermite 
polynomial Ha by 

2\a\ 



Ha{0 = -^e\^\'/'d^{e-\^\'/') , eeR^ 



where a\ = (ai!) ■ . . . ■ (a^l). It is not difficult to verify that Ha is a polynomial of degree 
|a| which satisfies C*Ha = — -^-ffo, where £* = — ■ is the formal adjoint of C. 
In addition, the following orthogonality relations hold for all a,f3 G N^: 

HaiOMO d^ = 5ap- 
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Finally, for any smooth function f{p), we have the relation 

(^^a(iv,)/) (0) = la; (/(p)el^l^) 



p=0 



Assume now that m > y. For all n G Z such that n + y < m, we define a continuous 
projection P„ : L^{m) L'^{m) by the formula 

|a|<n ^ ^ ' 

We also set Qn = 1 — Pn- Remark that P„ = and Qn = 1 if n < 0. If n > 0, it 
is clear from the definitions that Pn is the spectral projection onto the ^^=o 
dimensional subspace spanned by the eigenfunctions of C corresponding to the eigenvalues 
{— 1 1 A; = 0, 1, . . . , n}. For later use, we note that the condition P„/ = is equivalent to 

/ r/(0 = for all a e with |q;| < n . 

(4) The semigroup e^'~ . The operator C is the generator of a linear semigroup S{t) — e^'-' 
given by the following expressions: 

{SiP)l){p) = e-««l^lV^(pe--/2) , or 

(47ra(T))^ Jn^ ^ 4a(Tj / 

where a(r) = 1 — e~^. It follows from these formulas that S{t) is a strongly continuous 
(but not analytic) semigroup on L'^{m) for any m > 0. 

The main technical result of this section is the following estimate on the semigroup 
S{t) = e^^: 

Proposition A. 2 (a) Fix m > 0, and take n G Z such that n + y<m<n + l + y. 
For all a e and all e > 0, there exists C > such that 

\\d-S{r)QJ-\\m < ^(:^e5(f-+^)||/|U , (89) 

for all f e LP'im) and all t > 0. 

(b) Fix n e N U { — 1}, and take m e R such that m > n + 1 + For all a e and 
all e > 0; there exists C > such that 

||5"^(r)g„/|U < e-^ni/IU , (90) 

for all f e L'^{m) and all t > 0. 
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Remark A. 3 Estimate has been obtained in ^^j , j^, and /Q/ under the assumption 
that m is sufficiently large, depending on n. As is clear from Theorem \A.l\ , the condition 
m > n + 1 + Y is optimal. 

Using Proposition |A.2| , it is easy to complete the proof of Theorem |A.1| . Indeed, if 
m and n are as in part (a) of Proposition [A.2| , then (j{C) = (y{CPn) U a{CQn)- By 
construction, a{CPn) = if n < and a{CPn) = {0; — |; . . . ; — |} if n G N. On the other 
hand, by the Hille-Yosida theorem (see for instance |^), the bound ( |5UD imphes that 
a{CQn) C {A e C 1 3fJ(A) < f - f }. Thus 



a(£)c{AGc|5R(A)<^-|}u{ 



U 



and the reverse inclusion has already been established (see (1) and (2) above). Finally, 



since a{CPn) H ai^CQr. 



we see that the multiplicity of the eigenvalue A^ 



(fc = 0. 



is exactly { ) • This concludes the proof of Theorem |A.l 



QED 



Proof of Proposition |A.2| . We first show that (a) ^ (b). Let n G N U { — 1} and 

m > n + 1 + y. Taking fi G N such that n+Y<m<n+l + y, we can write, for any 
/ e L2(m), 



QJ = (1 - Pn)f = {Pn - Pn)f + Qnf . 



Since m — > n> n + 1, we have by 
WS{T)Q^f\U < 
On the other hand, since S{T)(j)c 



C 



e ~' 



a(r)l°l/2 



n<\a\<n 

ri + l 

we easily obtain \\d"S{T){Pn — Pn)f\\m < Ce"^""^ 
To prove (a), we first remark that, for all a G N^, 



r > . 



1/2 



(p^-p„)/= E (X^^-^^)^^^)^^ 

. This proves (|9 



(5°^(r)/)(0 



Nt 

e 2 



a(r) 



i(r)i/2 



/(e'e5)de' 



(91) 



where 0„ is defined in (|8|). Since 1 + l^]"^ < C(l + |^ - f |™)(1 + |f |™), we have 



;i + ier)|(9"5(r)/)(0l < c 



JVr 

e 2 



a r 2 



a(r)Vi 



Observing that 



a(r) 2 jRAf 



a r 



)l/2 



x(i + ie'r)i/(e'e5)idr . 

dr/ < C , 
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and applying Young's inequality \g * h\2 < \g\1\h\2, we obtain 

Nt 

WS{r)fU<C^^\\fU, r>0. (92) 
a(r) 2 

In view of (p2|), it is sufficient to prove ( pQ]) for r > 1. Since 1 — e~^ < a(r) < 1 when 
r > 1, we may drop for simplicity the factors a(T) in (|9T|). Thus, all we really need to 
prove is: 

Lemma A. 4 Let cj) G 5(R^), m>0, and n E Z such that n + ^ <m<n+l + Y- For 
any f G L'^{m), define 



(5(r)/)(0 



/ 0(^-^)/(^?e2)dr/ , r>0. 
J UN 

Then, for all e > 0, there exists C > such that, for all f G L'^{m), 

\\S{t)QJU < Ce-i(™+f-)||/|U , r > 0. (93) 
Proof: Expanding in Taylor series, we obtain for all C,,f] & R^: 

0(e-r])= 5^ ^^(9»(e)r^"+ J2 -^M^^^K. 

|a|<n |«|=n+l 

where ^ 

^a{^,v) = {n + l) [ (l-s)"(9»(e-sr/)ds. 

If = Qnf, then PnQ = 0, hence J-^n V°'9{v) = ^ with |«| < n. 

Therefore 



(^(r 



= (^i(r)^7)(0 + (^2(r)^7)(0 + (^3(r)^7)(0, 



where 



{Si{r)gm = / 0(e-r7)^7(r7e2)dr7 
'|r)|>i 

.1)1"! 



(^3(r)^)(0 = E <f«(e,r/)r^"^7(r/e^)dr/. 



|o|=n+l 



Since 1 + 1^1™ < <^(1 + 1^ - ^/Dl^/I" when |r/| > 1, we have 



(1 + \T)\{S^{r)gm\ <C (1 + le - - ^)ll^ri^7(^e^)l dv . 

'l'?l>i 
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Applying Young's inequality, we conclude that 

||(^i(r)(7)|U<Ce-5(™+f)||(7|U, r>0. (94) 
On the other hand, if |a| < n, we have by Holder's inequality 

ri''g{riei)dri < [ |r?r"'"(|r/r|^(r/e5)|) dr/ 

|»,|>1 i|r,|>l 

< ([ |r/|2("-™)dr/ ) e-5("+^)||^|U , 

Vi|r;|>l / 

hence 

||(^2(r)(7)|U<Ce-5("^+f)||(7|U, r>0. (95) 

To bound the last term S3{T)g, let a G such that |a| = n + 1, and define \I/q(^,?7) = 
(1 + \^\"^)^a{^,v)- r > 1, it is easy to verify that 

sup ( [ \^4^,r]Wdv) ' +snp([ \^^{^,r]WdA ^ < oo . 

As is well-known, this implies that the linear operator defined by the integral kernel 
\l/„(^,?7) is bounded from L^({|?7| < 1}) into Lp(R^) for all q < p. We now distinguish 
two cases: 

Case 1. If m < n + 1 (which is possible only if = 1), we have 



l + m\iS3ir)gm\ < J2 -J \^.{^,v)\\vr'\9ive^)\dv 

, , , , J\v\<i 



(96) 



\a\=n+l 



^ J2 -J \^a{^,v)\\vr\9{ve^dr^ ■ (97) 

Using (j9^) and the remark above with p = q = 2, we immediately obtain 

||(^3(r)^7)|U<Ce"i(-+f)||(7|U, r>0. (98) 

Case 2. If n+1 <m< ra+l+f , we define G [1, 2) by the relation A^(^-|) = m-{n+l). 
Given any e > 0, we choose q G (q'*,2) such that A^(^ — ■^) < e- As is easily verified, 
there exists C > such that ||'7|"^^fi'('7)|g < ||5'||m for all g G L'^im). Using ( P^D and the 
remark above with p = 2, we obtain 

\\{Hr)g)\\^ < c([ i\vr'\g{ve'^)\r dvY" 
\J\v\<i / 

< Ce~^("+^+f)||r^|"+i(^(r/)|, <Ce-^('"+f-^)||(7|U (99) 

since n + 1 + f = m + f - N{j^ - ^) > m + f - e. Combining (|9|), (0), and (H) or 
(PSp, we obtain (0). This concludes the proof of Lemma |A.4| , hence of Proposition |A.2 . 
QED 

The results of this section can easily be generalized to weighted spaces with p 2. 
For instance, the following bounds are used in Section |^ to control the nonlinearity in 
(0)- 



36 



Proposition A. 5 Let 1 < q < p < oo, m > and T > 0. For all a G N^, there exists 
C > such that 

Kd"S{r)f\, < — ^^--^ KfU , (100) 
a(r) 2 ^1 P>+ 2 

where b{0 = {1 + 

Proof: This estimate follows easily from (|9T|) and Young's inequality. QED 

B Bounds on the velocity field 

In this section, we study in more detail the relationship between the vorticity w{^) and 
the velocity field v(^) obtained from w(^) via the Biot-Savart law (|TB|). In particular, we 
obtain sharp estimates for the spatial decay of the velocity as |^| — > oo. This information 
is systematically used in Section § to derive properties of the solutions of the Navier-Stokes 
equation from the results we have on the vorticity. 

First, we give explicit formulas for the velocity fields corresponding to the first eigen- 
functions of the linear operator C defined in (|15|). For our purposes in Section ^, it is 
sufficient to consider only the first three eigenvalues Aq = 0, Ai = —1/2, and A2 = —1. 

1) The first eigenvalue Aq = is simple, with Gaussian eigenfunction 
The corresponding velocity field v'-^(,^) is given by 



1 e-l€lV4 -1/^2 



Remark that |v^(0| ~ 1^1"^ as |^| ^ 00, so that G L'^iR^Y for all q > 2. 

2) The second eigenvalue Ai = —1/2 is double, with eigenfunctions 

i^.(0 = 5.^(0 = -|g(0, ^ = 1,2. 

The corresponding velocity fields are 

v^HO = 5iv^(0 , v^^(0 = 92V^(0 . (102) 

Note that |v^'(0| ~ 1^1"^ as |^| 00, so that v^- G Li(R'^f for all g > 1. 

3) The third eigenvalue Ai = — 1 has multiplicity three. A convenient basis of eigenfunc- 
tions is 

H,{0 = AG(0 = (iieP-l)G(0, 
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and the corresponding velocity fields read 

1 



.Hi 



(0 = Av^(0 = 2G(0 



6 
-6 



(103) 



Remark that v^^ decreases rapidly as |^| oo. In constrast, for j = 2, 3, |v^^ (^)| ~ |^|^^ 
as lel ^ oo, so that v^^ G L^{R^y, but 6v^^ ^ ^^(R^)^. 

Next, assume that w G L'^{m) for some m > 3. Then w can be decomposed as follows: 

2 3 

w{0 = "^(0 + 5^ + ^7,i^i(0 + ^(0 , (104) 

i=i j=i 

where a, (3i are given by ( ^2]) and 7j by (|51|). The corresponding velocity field has a similar 
decomposition 

2 3 

v(0 = «v^(0 + E^^^'''(^) + E^^-^''^(^)+^''(0 , (105) 



where is obtained from R via the Biot-Savart law ([T6|). By construction, the remainder 
term i? in ( M satisfies /j^^ ^r^2'^(0 = for all ni, G N such that rii + na < 2. 
Using this information, we shall show that |v'^(^)| ~ 1^1"™" as |^| oo. In view of ( |104| ), 
it will follow that v G L^(R^)^ if and only if a = 0, namely 

f «;(Ode = 0. (106) 

Moreover, if (|106D holds, then v G L^{R?Y if and only if (3i = (32 = 0, namely 



R2 



iMO de 



6n;(0 de = 



R2 



(107) 

Finally, if (|106|) and ( p^ ) hold, then 6v G /.^(R^)^ if and only if 72 = 73 = 0, namely 

I i^!-^lMOd^= [ ei6^(Ode = 0. (108) 

JR2 JR2 

Remark that it is not necessary to assume here that 71 = 0, since the velocity field v^^ 
(unlike v'^2 a^^j v^^) decreases rapidly at infinity. 

Proposition B.l Let w G L'^{m) for some m > 0, and denote by v the velocity field 
obtained from w via the Biot-Savart law f|7^. Assume that either 

1) < m < 1, or 

2) 1 < m < 2 and ( \10(^ holds, or 

3) 2<m<3 and ( \iM ), (MD ^old, or 

4) 3 < m < 4 and (^UB^, (^U^, ( \TU8i ) hold. 

If m ^N, then for all q G (2, +00), there exists C > such that 



(109) 



where b{^) 
then 



1 + l^n^^^- IfmeN and b'^w G ^^(R^) n ^'^(R^) for some p < 2, r > 2, 



(110) 
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The proof of Proposition B.l relies on the following weighted Hardy-Littlewood- 
Sobolev inequality: 

Lemma B.2 If < m < 1 and 

then, for all q E (2, +oo), < C\b"'uj\2. 

Proof: We use the dyadic decomposition 

oo 
j=0 

where Bq = e \ \^\ < 1} and Bj = g | 2^'^ < \^\ < 2^} for j e N*. Let 
Ui = ulsi and = wls^, z G N. Clearly Ui = J^jeN ' where 



Jb, Is ~ V\ 



Fix q E (2, +oo), and define p E (1, 2) by the relation ^ = ^ ~ |- 

If K j| — 1) follows from (|D that |A.y|g < Clcu^lp. By Holder's inequality, \ujj\p < 
mes{Bjy/''\ujj\2 < C2^^/'^\ujj\2. Thus, |A,j|, < C2^^/^\ujj\2. 

If K ~il — 2, Young's inequality implies that \^ij\q < M]^~'''^'^M|'''^|(X'j I2, where where 
Mi = sup / -dr/ , M2=sup / -de . 

If 2 > j + 2, then 1^ - r/| > |e| - jr^] > 2*"^ - 2-'' > 2^'^ for all ^ E Bi, r] E Bj. Thus 
Ml < C2-^mes(5j)^/P < C2-'2'^^/p and M2 < C2-*mes(Si)^/P < ^2-^22^/^, for some 
C > independent of ij. It follows that \Aij\g < C2''^'-^h^'/'i\ujj\2. If j > i + 2, then 
le — r/l > 2^~^ for all ^ G -Bi, G -Bj, and a similar calculation shows that \Aij\q < 
C2'^^^'^\ujj\2- Summarizing, we have shown that 

2-^'/''\u,\q<Cj2Kij\u;j\2 , tEN, 
jeN 

where Kij = 2~5l*~-?'l~5(*~-?). Now, by definition of the sets Bi, there exists C > 1 such 
that C-12™ < 6(0 < C2™ for all ^ E Bi and all i G N. It follows that 

where ir^""^ = 2-3N-il+('»-5)(*-i). Since < m < 1, < 2"°l^"Jl for some a > 0, 

hence K^"^^ defines a bounded linear operator from £^(N) into i''(N). This concludes the 
proof. QED 
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Proof of Proposition B.l . The proof is naturally divided into four steps, according to 
the values of m. 

1) If V and w are related via the Biot-Savart law (|16|), then 



|v(OI < 



2tt 



1 



R2 



\w{T])\dr], ^gR2 



[111) 



and (|109|) follows immediately from Lemma [B.2| . To prove ( [1101 ), we remark that b{C,) < 
1 + 1^1 < 1 + 1^ - ^1 + |^?| < 1^ - + 26(7]) for all ^,r] e hence 



KO|v(OI < 



27r 



|w(?7)| di] + 



R2 



R2 



■b{ri)w{ri) drj 



(112) 



In view of (|^), the second integral in ( [L12| ) is uniformly bounded by C{\bw\p + \bw\r) if 
p < 2 and r > 2. On the other hand, it is clear that \w\i < C\bw\p. This concludes the 
proof of case 1. 

2) For all ^, 77 G with ^ 7^ and ^ 7^ r], we have the identity 



^1 - Vi 6 



(ei-^i)(|-r/) + (6-^2)(^ Ar7) 



where ^ ■ ^7 = ^1771 + ^2^2 and ^ Arj = ^irj2 — ^2^1- Therefore, for all w G L'^ijn) satisfying 



TO, we find 



1^2(01 = ^ 

ZTT 
< C 



R2 



- ^ ) w{ti) dr] 



1 



R2 



rj\\w{ri) \ dr) . 



A similar bound holds for the first component fi(0 of the velocity field. Combining these 
estimates with ( p. 1 1| ) , we obtain 



m^m<c 



R2 



■\b{r])w{T])\d7] , ieYC . 



This inequality has exactly the same form as (|111|) , with v, w replaced by bv,bw. There- 
fore, applying the preceding result to 6v, bw instead of v, w, we obtain ( |109| ) and ( |110[ ) in 
case 2. 

3) For all ^, G R^ with 7^ and ^ 7^ r/, we have the identity 



^—j, ((^1 - Vi)m ■ rif - mv?) + 2(6 - V2){i ■ r7)(C A ry)) . 



Therefore, for all w G L'^{m) satisfying ( |106|) and ( |107] ), we find 



1^2(01 = TT 

ZTT 

< c 



R2 



1 



^(77) d77 



R2 \m-V\ 



7]] \w{ri)\ df] . 
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A similar bound holds for the first component Vi{^), hence 

1 



m)Mo\<c 



R2 



■\b{r]yw{r])\dr] , ^ G 



This inequality has again the same form as (|111| ), with v, w replaced by 6^v, b'^w. There- 
fore, proceeding as above, we obtain (|109|) and ( |11CI| ) in case 3. 

4) For all ^, ?7 G with ^ 7^ and ^ ^ rj, we have the identity 



6 -Vi 6 - Vi 



\C-v\ 



1 



(113) 



+ {^2-V2){^Arjm^-rjf-\^\'\r]\ 



The left-hand side of ( |113|) has the form {^i—Tii)/\$,—ri\'^ + A{^, t]), where A is a polynomial 
of degree 2 in 77 with ^-dependent coefficients. As is easily verified, the terms in A that 
are quadratic in rj can be written as 

For all w G L'^{m) satisfying ( |106| ), (|107|) and ( |108|) , we thus have Jj^a A(^, ri)w{ri) di] = 0. 
Therefore, using (|113|) , we obtain 



?7|^|w(r7)| di] . 



A similar bound holds for the first component vi{C,), hence 

1 



m)MO\<c 



R2 



\b{riYw{ri)\dri , ^ G R^ . 



This inequality has again the same form as ( |111| ), with v, w replaced by 6^v, b^w. There- 
fore, proceeding as above, we obtain (|109|) and ( |11CI| ) in case 4. This concludes the proof 
of Proposition |Kl[ QED 

Corollary B.3 Assume that w G L'^{m) for some m > 1, and denote by v the velocity 
field obtained from w via the Biot-Savart law I^Td{). Then v G L^(R^)^ if and only if (\1U(\) 
holds. 

Proof: Without loss of generality, we assume that 1 < m < 2. For any w G L'^{m), 
we have the decompositions w = aG + w, \ = av*-^ + v, where a = ""^(O and w 
satisfies the assumptions of point 2 in Proposition [B.l| . Setting q = 2/(m— 1) in ( |109| ), 
we obtain 6v G L'^(R^)^, which implies v G L^(R^)^. Since v*^ ^ L^(R^)^, it is clear that 
V G L2(R2)2 if and only if a = 0. QED 

Corollary B.4 Assume that w G L'^{m) for some m > 2, and denote by v the velocity 
field obtained from w via the Biot-Savart law ^Td{). Then v G ^-'^(R^)^ if and only if (\10(^ 
and l\10V hold. In that case, d^ = for i = 1,2. 
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Proof: Without loss of generality, we assume that 2 < m < 3. For any w G L'^{m), we 
have the decompositions 



w = aG + (3iFi + (32F2 + w , V = av^ + /^iv^^ + (3^^^^ + v 



where a,/5i,/32 are given by (|4^) and to satisfies the assumptions of point 3 in Proposi- 
tion Thus fo^v e with g = 2/(m-2), hence v e L^(R2)2. On the other 



hand, it is easy to verify that v — v G L^{R?Y only if a = /?i = = 0. Thus, 

V G L^(R^)^ if and only if ( p.06| ) and ( |107] ) hold. In that case, since divv = 0, we must 
have Vi{i) d^ = for i = 1, 2. QED 

Corollary B.5 Assume that w G L'^{m) for some m > 3, and denote by v the velocity 
field obtained from w via the Biot-Savart law (|7^. Then h\ G L^i^Y ^/ ^'^^ ^^^V ^/ 
i^Md[ ), ^TOlj ) and ( ^TOBj ) hold. In that case, 



[ ei^i(Ode = / 6^2(0 = 0, 

/ 6t^i(0de = - / 6^2(0 de = 71 • 

7r2 Jr2 



Proof: Without loss of generality, we assume that 3 < m < 4. Let w G L'^{m), and 
consider the decomposition ( |104| ), ( |105| ). Then R satisfies the assumptions of point 4 in 
Proposition |B.1| . In particular, 6^v^ G L''(R^)^ for q = 2/(m— 3), hence G L^(R^)^. 
On the other hand, it is easy to verify that 6(v — v^) G L^(R^)^ if and only if a = /3i = 
/52 = 72 = 73 = 0. Thus, 6v G Li(R2)2 if and only if (|106D, ([lO^) and (|108D hold. 

Assume now that 6v G L^(R^), and consider the vector field A(^) = .^^v(^). Then 
A G L2(R2)2 and div A = G ^^(R^)^, hence 



0= / divA(Ode= / 2^M0di. 

7r2 Jr2 



Similarly, using the vector fields and (,1^2^, one obtains ^2^2 d^ = and J-^2{^iV2 + 
^2Vi)d^ = 0. Finally, if B(0 = |^pv^ = \^W~V2,Vi), then B G L'^iR^ and divB = 
-l^l'^w + 2(^2^1 - ^iV2) G ivi(R2)2, hence /^^^ divBd^ = 0. Using (^, we thus find 

471=/ |er^(0d^ = 2/ (6^i(0-ei^2(0)de = 4 / 6^1 (Od^ 

JR2 JR2 JR2 



This concludes the proof of Corollary |B.5| . QED 



C Proof of Proposition |4.10| 

The existence and attractivity of the manifold follow from the results of exactly as in 
the proof of Theorem |3.5| and Corollary By construction, the function Q whose graph 
gives the invariant manifold satisfies the integral equation (see, for instance, [p]] ): 

/■oo 

gif3, 7) = / e-^AAo(e-/2/3, 0-^(7; g{e^^'(3, 0-^(7; /3))) dr , (114) 
^0 
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where (/3, 7) 1— > (e"""/^/?, 0^(7; /3)) is (the projection onto of) the flow defined by 
equations ( |60|) on Wjf"^, and A/q is the nonhnear term in the equation for p in (|60|) , cutoff 
outside a neighborhood of the origin as in p3|). Note that in computing 0'^(7; f3), one can 
first solve the equations for f3 explicitly and then insert these solutions into the equations 
for 7j, thus treating them as inhomogeneous terms. 

Remark C.l Note that the nonlinear terms fj in the equations for 7j are also cutoff so 
that they vanish outside a neighborhood of the origin of size tq. It is an easy exercise to 
show that for any e > 0, there exists tq > sufficiently small so that the resulting flow 
4>'^{'j; f3) satisfies 

|0-^(7;/3)|<Ce(^+^)-(|7| + |/3p) , (115) 

for all T > 0. 
Define a mapping 

POO 

H'3mi) = / e-^Aro(e-/2^,r"(7;/3),6;(e^/2Ar"(7;/3)))dr . (116) 

Equation ( |114| ) implies that has a fixed point in the space of functions with globally 
bounded Lipshitz constant. Moreover, if one fixes the cutoff function Xro in (|23|) , then 
(for ro > sufficiently small) is a contraction in this function space, so that Q is the 
unique fixed point of J^. Thus, the proof of the proposition is completed if we can show 
that for some C5 > 0, ^ maps the set of functions satisfying (^) into itself. 

Let X be the Banach space x x Eg, equipped with the norm || (/3, 7, p) ||x = 
+ I7I + IIpIU' let Xt-o be a cutoff function like that introduced just prior to (pSf) , 
but defined with respect to the norm || ■ ||x rather than || ■ Then the nonlinear term 
in (|114|) can be written as 

Aro(/3, 7, p) = a{V ■ (x.o(/3, 7, p)F(/3, 7, p)) - • V(xro (/?, 7, P)p)} , (117) 

where Qs is the projection onto Es- Commuting the derivatives in the nonlinear term 
through the factor of as we did in (|19D, the right-hand side of (|116|) can be rewritten 

as 

'aV-e-(^-^)[(x.oF)(/3(r),7(r),^(/3(r),7(r))) 



_ve(/3M,.M)^^^^(^(^)^ ^(/3(r), 7(r)))^(^(r), 7(r))] dr 

/oo 
e(--i)^Q,V ■ e^-^ [(x.oF)(/3(r), 7(r), ^(/3(r), 7(r))) 

_^gmr)Mr)) ^^^^^(3^r),j{T), ^(/3(r), 7(r)))^(/3(r), 7(r))] dr , 

where in keeping with our usual notational convention, denotes the velocity field 
associated with the vorticity field Q. Also, to save space, we have used the short-hand 
notation /3(t) for and 7(t) for ^""^(7; (3). Fixing e > 0, we can take the || ■ II4 norm 

of this expression and use the estimates on the semigroup from Appendix ^ to obtain 

/■oo 

\\H<3mi)\U <C e-^i-^)(^-i)[(x.o||F||4)(/3(r),7(r),^(/3(r),7(r)))||4 

+ x.o(/3M, 7(r), ^(/3(r), 7(r))) ||^(/3(r), 7(r)) dr , (118) 
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where we estimated the term involving v^(xro(/3, 7, 7))^) in the same way as we 
estimated the term involving v(s)w(s) in Lemma Using the estimate on F in (|57D , 
we see that 

||Xn,(/3,7,^)F(A7,^)l|4 < Cm\ + I7l)(l7l + II^IU) + ||^||Dx.o(A 7, ^) • 

If we now assume that ||^(/3,7)||4 < Cs{\P\^^^ + 171^^*^), then, using the fact that 
Xroif^yly p) = if either \/3\ or I7I is bigger than r^, it is easy to verify that 

\\Xroif3, 7, G)F{P, 7, Q)\U< C4{\(3t' + . (119) 

In a similar fashion we can bound 

xMi,Gif3n))mf3,iml < c4i\f3t' + l7l^~') ■ 

Inserting these estimates into (|1 18|) and using (|115| ) to bound 7(r) = ^^^(7; /3), we find 

POO 

mGm-f)\U < Crd e-(i-^)(--i)[e(3-^)(-/2)|/3|3-^ + e(i-^)(i+^)-|7|^-']dr 

Jo 

< Cie,6y,i\f3t' + h\l-') , (120) 

provided that e < 25/5. For fixed e and 6 (satisfying e < 25/5) we can always pick tq 
sufficiently small so that C(e, 6)rQ < Cs- Thus JF maps the set of functions satisfying (p2| ) 
into itself (for ro sufficiently small) and hence the fixed point of J-' whose graph gives the 
invariant manifold must satisfy (|62D. This completes the proof of Proposition |4.1CI| . QED 
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